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Abstract 

We report the first determination of the relative strong-phase difference between D° — > 
Ks L K + K~ and D° — > K$ L K + K~ . In addition, we present updated measurements of the relative 
strong-phase difference between D° — > Kg ' and D -»• K° s r iT + iT . Both measurements ex- 

ploit the quantum coherence between a pair of D° and D° mesons produced from ^(3770) decays. 
The strong-phase differences measured are important for determining the Cabibbo-Kobayashi- 
Maskawa angle 7/^3 in B~ — > K~ D° decays, where D° is a D° or D° meson decaying to Kgh + h~ 
{h = 7T, K), in a manner independent of the model assumed to describe the D° — > Kgh + h~ de- 
cay. Using our results, the uncertainty in 7/^3 due to the error on the strong-phase difference is 
expected to be between 1.7° and 3.9° for an analysis using B~ — > K~D°, D° — > K^n + ir~ decays, 
and between 3.2° and 3.9° for an analysis based on B~ — > K~D°, D° — > KgK + K~ decays. A 
measurement is also presented of the CP-odd fraction, of the decay D° — > K^K + K~ in the 
region of the <p — > K + K~ resonance. We find that in a region within 0.01 GeV 2 /c 4 of the nominal 
<j) mass squared T- > 0.91 at the 90% confidence level. 



*Now at: Pacific Northwest National Laboratory, Richland, WA 99352 
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I. INTRODUCTION 

A central goal of flavor physics is the determination of all elements of the CKM matrix pQ , 
both magnitudes and phases. Of the three angles of the b — d CKM triangle, denoted a, (3, 
and 7 by some, (fi 2 , <fii, and 3 by others, the least well determined is 7/03, the phase of 
relative to V u b. It is of great interest to determine 7/^3 using the decay B ± — y K ± D°, since 
in this mode, the 7/(^3 value is expected to be insensitive to new physics effects in B decay; 
here, D° is either a D° or D meson decaying to the same final state. This is in contrast with 
most measurements of CP violation, which are dominated by processes that have significant 
contributions from loop diagrams that can be influenced by new physics [21 EJJ. Therefore, 
precise measurements of 7/^3 from the decay B ± — y K ± D° compared to the predictions for 
7/03 from loop-dominated processes provide a stringent test of the origin of CP violation 
in the Standard Model. Sensitivity to the angle 7/03 comes from the interference between 
two Cabibbo-suppressed diagrams: b — y cus, giving rise to B~ — y K^D , 1 and the color and 
CKM suppressed process b —y ucs, giving rise to B~ —y K~D jlj. Promising D° decays 
for measuring 7/(^3 using this method are D° —y KgTr + 7r~ E] and D° —y KgK + K~ , 
here designated collectively as D° —y Kgh + h~ . To make use of these decays, however, it is 
necessary to understand the interference effects between D° —y K®h + h~ and D° —y K° s h + h~ . 
These effects can be measured using CLEO-c data. A study of the decay D° — y Kgn + n^ 
has already been published [7] . Here we present an update of that analysis, and first results 
from the decay D° -> K° S K + K~. 

Let us write the amplitude for the B~ —y K~D°, D° —y K^h + h~ decay as follows: 

f B -(m 2 + ,m 2 _) ex f D (ml,m 2 _)+r B e^ B -^fn(ml,m 2 _). (1) 

Here, m 2 + and m 2 _ are the invariant-mass squared of the K^h + and Kgh~ pairs from the 
D° — y Kgh + h~ decay, which define the Dalitz plot, rn 2 _)(fj)(m 2 l _, m 2 _)) is the am- 

plitude for the D°(D ) decay to K° s h + h~ at (m 2 ] _,m 2 _) in the Dalitz plot, tb is the ratio 
of the suppressed to favored amplitudes, and 5b is the strong-phase difference between the 
color-favored and color-suppressed amplitudes. Ignoring the second-order effects of charm 
mixing and CP violation in charm [5j[8j|9], we have fj){m 2 + , m 2 _) = fc(m 2 _ } m+), and Eq. (JTJ 
can then be written as: 

f B -(m 2 + ,m 2 _) ex f D (m 2 + ,m 2 _)+r B e^ B -^f D (m 2 _,m 2 + ). (2) 

The square of the amplitude clearly depends on the phase difference A5d = SD(rn 2 h ,m 2 _) — 
Sd^ 2 -, where 5i)(m^_,m?_) is the phase of /d(j71^_, m^.). Thus, for the determination 
of 7/03, one must know A5d- 

Analyses of B~ —y K~D° decays to date extracted A5c>(m;j_, m 2 _) for each final state by 
fitting flavor-tagged D° ^K° s n + n- [TUHH] and D° -^K S K + K~ ^ E] Dalitz plots to D°- 
decay models involving various two-body intermediate states. The systematic uncertainty 
associated with this modeling is hard to estimate; the assigned values vary between 3° and 
9° for the more recent analyses. In order to exploit fully the high statistics expected at 
LHCb [T5l [16] and future e + e~ 5-factory experiments [HI [18] it is highly desirable to avoid 



Here and throughout this paper the charge-conjugate state is implied unless otherwise stated. 
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this modeling uncertainty, and to do it in a manner which keeps all other error sources small 
compared with the foreseen statistical precision. 

In the analysis presented here, we employ a model-independent approach to obtain 
A#D(m^, m 2 _) as suggested in Refs. [19], by exploiting the quantum coherence of D° — D 
pairs at the ?/>(3770). Because of this quantum correlation, K$h + h~ and K < [h + h~ decays re- 
coiling against flavor tags, CP-tags, and D° — > Kgh + h~ tags, taken together, provide direct 
sensitivity to the quantities cosA^d and sinA^ for each final state. The analysis is per- 
formed in discrete bins of D° —>Kgh + h~ Dalitz space. We have updated the D° -^■Kg7i + 7i^ 
analysis reported in Ref. [7] by providing results in alternative sets of Dalitz-plot bins, and 
by reducing some of the systematic uncertainties. 

In addition measurements of the time-dependent evolution of the D° — > Kgh + h~ 
Dalitz plot provide some of the most precise constraints on charm-mixing parameters [20J. 
These measurements also rely on P°-decay models that introduce significant systematic 
uncertainties. A model-independent determination of the charm-mixing parameters from 
D° — >■ Kgh + h~ that uses the same strong-phase difference parameters as the 7/^3 analysis 
of B~ decay has been proposed [9]. The advantage of the model-independent approach is 
again the elimination of model-dependent assumptions about the strong-phase differences. 

We also present the first model-independent measurement of the CP content of the decay 
D° — > KgK + K~ in the region of the — > K + K~ resonance. The decay D° — > K®(f) is a CP- 
odd eigenstate and has been used as such in several analyses; see for example Refs. [2U [22] . 
The cf) — > K + K~ resonance is usually defined by a mass window about the nominal mass. 
Despite its narrow natural width of 4.26 MeV/c 2 [23], the potential contributions from 
CP-even final states beneath the (f> resonance, such as D° — > Agao(980) and non-resonant 
D° — > KgK + K~ decays, must be accounted for. Using D° — > Kg L K + K~ decays recoiling 
against CP eigenstates we determine the CP-odd fraction of decays, in the region close 
to the (f) resonance. A measurement of J-L allows a systematic uncertainty related to the 
CP-even contamination to D° — > K®<j) decays to be assigned without assuming an amplitude 
model for the decay D° -> K° S K + K~ . 

This paper is organized as follows. The formalism for the measurement of the strong-phase 
difference and T~ is outlined in Sec. [TTJ The choice of Dalitz-plot bins is given in Sec. |III[ 



The event selection is described in Sec. [TV] Sections [V] and [VI] present the extraction of 
the variables associated with the strong-phase differences and the assignment of systematic 
uncertainties, respectively. The impact of these results on the measurement of 7/^3 is 



discussed in Sec. VII, along with the measurement of A summary is given in Sec. VIII 
Throughout this article the D° -^■Kgir + 7r^ and D° —^KgK + K~ analyses are described in 
parallel, but more weight is given to the latter as it has not been presented previously. 



II. FORMALISM 

Giri et al. proposed (3] a model-independent procedure for obtaining A5^(m^_,m 2 ) as 
follows. The Dalitz plot is divided into 2J\f bins, symmetrically about the line m 2 + = m 2 _. 
The bins are indexed with i, running from —M to M excluding zero. Thus the coordinate 
exchange m 2 + <H- m 2 _ corresponds to the exchange of the bins i <H- —i. The number of events 
(Ki) in the i th bin of a flavor-tagged Kgh + h~ Dalitz plot from a D° decay is then expressed 
as 

Ki = A D [\f D {ml,m 2 _)\ 2 dm 2 + dm 2 _ = A D Fi, (3) 
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where the integral is performed over the i th bin. Here Ad is a normalization factor and Pj 
is the fraction of D° — > Kgh + h~ events in the i th bin. The interference between the D° and 
D° amplitudes is parameterized by two quantities 



^ = / \f D (m + ,m_)\\f D (m_,m + )\ cos[A8 D (m + ,m_)]dm + dm_, (4) 

and ^ 

si = / \f D (m\,m^)\\f D (mt,m 2 + )\ sm[A5 D (m+,mt)]dm\dm 2 _. (5) 

\/fit_i Ji 

The parameters q and Sj are the amplitude-weighted averages of cos A5d and sinA^u over 
each Dalitz-plot bin. 

Though the original idea of Giri et al. was to divide the Dalitz plot into square bins, 
Bondar and Poluektov noted [19] that alternative bin definitions will lead to significantly 
increased sensitivity. In particular, one can choose to minimize the variation in A5e> over each 
bin according to the predictions of one of the models developed on flavor-tagged data [TUl - 
IT4"] . Note that this approach does not introduce a model-dependence in the final result for 
7/^3. This result will remain unbiased by the choice of an incorrect model, but will have 
less statistical sensitivity than expected. If we divide the Dalitz plot into M bins of equal 
size with respect to A5d as predicted by one of these models, then in the half of the Dalitz 
plot m 2 + < m 2 _, the i th bin is defined by the condition 

2n(i - 3/2) /j\f < A5 D (m 2 + , m 2 _) < 2ix{% - 1/2) /j\f, (6) 

and the — i th bin is defined symmetrically in the lower portion of the Dalitz plot. The choice 
of D° -^■KgTT + 7T~ binning with Af = 8 as obtained from the model presented in Ref. [12] is 
shown in Fig. [TJ A discussion on alternative choices of binning for D° — >Kg7r + 7r~ and those 



for D° ^K° S K + K- can be found in Sec. |TTT 



We now describe how CLEO-c data can be used to determine q and Sj. The event 
yields in the i th bin of both flavor-tagged and CP-tagged D° — > Kgh + h~ Dalitz plot are 
required. Because the ^(3770) has C — — 1, the CP eigenvalue of one D meson can be 
determined by reconstructing the companion D meson in a CP eigenstate. With a CP- 
tagged D° — > Kgh + h~ decay, the amplitude is given by 

fcp±(m 2 + ,m 2 _) = -^=[/ D (m+,m 2 ) ± f D (m 2 _,m 2 + )}, (7) 

for CP-even (+) and CP-odd (— ) states of a D° -^-Kgh + h~ decay. Since the event rate 
is proportional to the square of this amplitude, the number of events in the i th bin of a 
CP-even or CP-odd tagged Dalitz plot is then 

Mf = h CP± (Ki ± 2c^K l K. l + K_ i ), (8) 

where hcp± = S ± /2Sf is a normalization factor that depends on the number, Sf, of flavor- 
tagged signal decays, and the number, S^, of P°-mesons decaying to a CP eigenstate in the 
sample irrespective of the decay of the other D meson; this is referred to as a single-tagged 
(ST) sample. Alternatively the normalization factor can be defined in terms of branching 
fractions, £>, as h C p± = Bcp±/2Bf, where Bcp± (Bf) is the branching fraction of D° to CP 
eigenstates (flavor tags). Thus, access to q is enabled by measuring the number of events, 
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FIG. 1: Equal Afe binning of the L>° -^Kgvr+vr- Dalitz plot with J\f = 8 based on the model from 
Ref. |12j . The color scale represents the absolute value of the bin number, 



M? 1 , in a CP-tagged K^h + h Dalitz plot, and the number of events, in a flavor-tagged 
K%h + h- Dalitz plot. 

Important additional information can be gained through analysis of D° -^-Kgh + h~ vs. 

D° -^-Kgh + h" data. The amplitude for ^(3770) decaying to a double Kgh + h~ final state 
is as follows: 



2 2 /2 , 2 f D {m 2 + ,m 2 _)f D (m' 2 ,m'l) - f D (m%m' 2 )f D (m 2 _,m 2 + ) 
f(m + ,m_,m + ,m_) = ± ^- y= ± ^. 



(9) 



The primed and unprimed Dalitz-plot coordinates correspond to the Dalitz-plot variables of 
the two D° -^-Kgh + h~ decays. Defining My as the event rate in the i th bin of the first and 
the j th bin of the second D -^-K° s h + h~ Dalitz plots, respectively, we have: 

(10) 



Mi. 



Here, h corr 



N D jj/2Sj = N D jj/8B 2 , where N D j) is the number of DD pairs, and as before 
Sf is the number of flavor-tagged signal decays. Thus analysis of both D° -^-Kgh + h~ vs. 
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D° -^-Kgh + h~ data and CP-tagged D° -^-Kgh + h~ decays allows q and s, to be determined. 
The ambiguity in the sign of s^ can be resolved using weak model assumptions. 

The decay D° — >K^h + h~ , due to its close relations with D° —>Kgh + h~, can be used 
to improve further the q and Sj determination. We assume the convention that A(D° — > 

K° s h + h-) = A(D K° s h-h+). Then, since the Kg and K\ mesons are of opposite CP, it 
follows that A(D° K° L h + h-) = -A(D° Klh~h + ). Hence for K° L h + h~ the Dalitz-plot 



rates of Eqs. (|8j) and (10) become 

M't = h CP± {K> =F 2^K>K>_ l + K'_ l ) ) (11) 

and 



Mtj = hcorriKiKLj + K-iK'j + 2JK i K'_ j K_ i K' j (c i c , j + s^)), (12) 



for CP vs. D° -^K° L h + h~ and D° -^K s h + h~ vs. D -^-K Q L h + h~ respectively, where K[ c- 
and s\ are associated with the decay D° -^■K < [h + h~ . In the analysis all four parameters q, 
Si, c'j, and s[ are determined for each channel, but in order to improve the precision on q 
and Si constraints are imposed on the differences Aq = c£ — q and As, = s' { — Sj. These 
constraints are discussed in Sec. IVl 

Because the branching fraction of D° — > K°7T + 7r~ is around five times larger than D° — > 
K°K + K~ it is advantageous to first determine the coefficients q, Sj, c£, and for the 
former decay, and then use these to help improve our knowledge of the coefficients for 
D° K°K + K~ . This is achieved through measuring the bin-by-bin rates for K^K + K~ 
vs. KgTr + 7c~, K®K + K~ vs. KgTT + ii~, and K$K + K~ vs. K'jix+ix^ Dalitz plots, and using 



suitably modified forms of Eqs. (10) and (12). 



The expression for the CP-odd fraction in the region of the — > K + K resonance in 



D° KlK + K~ decays follows from Eq. (ph. We note that Mf + M~ = M+ + M. 



m 



addition, this sum is proportional to Ki + K_i, the total rate of D° — > K^K + K~ decays 
in the combined % and — % bins. Therefore, if bin % defines an interval of K + K~ invariant- 
mass squared, m 2 K+K -, about the nominal mass squared, then the CP-odd fraction of 
D° — > KgK + K~ decays, J--, in that region is 

j = Mr + MZl 

+ MZi + M+ + M+' K J 

We determine T- for four different invariant-mass squared intervals: 0.006, 0.010, 0.014, 
and 0.018 GeV 2 /c 4 . 



III. DALITZ PLOT BIN DEFINITIONS 

Measurements of cp and sj are presented for three and four alternative binnings for 
D° — > K$K + K~ and D° — > Kg7r + 7c~, respectively. The motivation for these choices and 
the resulting binning definitions are presented in this section. 



A. Binnings of the D° ->- K$K + K~ Dalitz plot 

We use an amplitude model determined by BABAR jH] for which a lookup table of 
the results in Dalitz space has been provided by the authors [21]. The amplitude model 
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uses the isobar formalism and consists of eight intermediate resonances, of which five are 
parameterized with Breit-Wigner lineshapes and three, ao(980)°Kg and ao(980) ± ii' =F , are 
parameterized by a coupled-channel Breit-Wigner function |25j . 

We consider binnings in which the Dalitz plane is divided into M = 2, M = 3, and 
M = 4 equal A8b bins, according to Eq. (J6]). A smaller number of bins provides superior 
statistical precision on the parameters associated with A5d but a reduced sensitivity to 7. 
Using a larger number of bins is not feasible due to the limited statistics available in the 
CLEO-c data; the fit to the parameters (Sec. [V]) fails to converge if J\f > 4. However, these 
alternatives will allow flexibility in matching an appropriate number of bins for the size of 
the available 5-decay sample when the values Cj and Sj are used to extract 7. 

The three cases considered are shown in Fig. [2] In each case, there is a narrow bin located 
at low values of m 2 K+K -, which is close to the diagonal boundary of the Dalitz plot. This 
bin encompasses the <ft intermediate resonance and typically contains the largest number 
of events. For three and four bins there are 'lobes' at high values of m 2 K+K - that contain 
relatively few events. 



B. Binnings of the D° 



Dalitz plot 



The four binnings used in the updated analysis of D° — > Kq7t + tt are described in this 



section. 



Sec. IIIB1 



The BABAR model [12] that is used to define the bin choices is described in 
Then the binning in equal intervals of the strong-phase difference is presented 



in Sec. Ill B 2 In Sec. Ill B 3 the procedure to optimize the binning for maximal sensitivity 



to 7 is described and the resulting binning is presented. In Sec. |IIIB4| we describe a modified 
procedure of optimization which takes into account expected background levels at LHCb. 
In Sec. Ill B 5 the binning in equal intervals of the strong-phase difference as defined by 

is given. Finally, in Sec. Ill B 6 we summarize the differences 



the latest Belle model 

between the various bin definitions and also assess what consequences a very recent BABAR 
D° —7- Kg7r + 7T~ model [14J, not available originally for our analysis, would have for the bin 
definitions. 



1. BABAR K-matrix model 

The amplitude models used by Belle PHI [32] and the first BABAR analysis [TT] are 
parameterized in terms of a Breit-Wigner isobar model. However, the broad tttt and Kit 
S'-wave components are not well described by such Breit-Wigner lineshapes. In particular an 
additional intermediate state, a', is required to fit the 7T7t 5*-wave even though it is known not 
to be a physical resonance. Furthermore, parameterizing these broad overlapping states in 
terms of Breit-Wigner functions is unphysical in that unitarity can be violated. Therefore, 
a more recent BABAR model [12] uses the .ff -matrix [2E] ansatz to parameterize the nir 
S-w&ve contributions, which does not violate unitarity. In addition, the empirical LASS 
lineshape [27] of the ^(1430) is used to improve the fit to the Kn S'-wave component. We 
refer to this approach as the BABAR 2008 model. The reduced x 2 for the BABAR 2008 
model fit to the D* + ->• D°ir + data is significantly improved over the first BABAR model 
[TT] and is much better than that for the Belle model. Therefore, this model was considered 
the best available and is used to define the nominal binnings implemented in this analysis. 
The results of the model have been made available by the BABAR collaboration in the 
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FIG. 2: Equal A5 D division of the D° -)• K^K+R- Dalitz plot into (a) J\f = 2, (b) M = 3, and 
(c) M = 4 bins. 



form of a lookup table of the amplitude and strong-phase difference, 5d, in a fine grid of 
0.0054 GeV 2 /c 4 x 0.0054 GeV 2 /c 4 sub-bins of the Dalitz plot variables m\ [21]. The m\ 
resolution, estimated from simulation, is of the same order as the sub-bin dimensions; the 
resolution is 0.006 GeV 2 /c 4 (0.015 GeV 2 /c 4 ) for D° -> K^tt+tt' (D° ->■ ^tt+tt-). (The 
lookup table for D° -> K° S K+K- has a grid of 0.0018 GeV 2 /c 4 x 0.0018 GeV 2 /c 4 sub-bins. 
The resolution is 0.005 GeV 2 /c 4 (0.010 GeV 2 /c 4 ) for D° -> K° S K + R- (D° ->■ K L K + R-).) 

Since performing the analysis of the CLEO-c data using the bin choices described here, 
which are based on the BABAR 2008 model, a new BABAR measurement [H] has been 
published which presents an updated version of the D° — > K®it + tt~ decay model that we 
term the BABAR 2010 model. This model is derived from a larger data sample that has 
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been reprocessed and the analysis has been improved with respect to experimental systematic 
uncertainties. In Sec. Ill B 6 we assess the consequences on the bin choices of the differences 
between the two models and conclude that they are minor. 



2. Equal A5d binning of the BABAR 2008 model 

The first binning of the Dalitz space for the BABAR 2008 model we consider follows the 
proposal in Ref . [19] , which was used in the previous CLEO-c analysis [7] and in the analysis 
of D° — > K$K + K~ . The binning is such that there are M = 8 bins of A5d in each half the 
Dalitz plot as defined in Eq. ([6]). This equal A5d binning for the BABAR 2008 model is 
shown in Fig. [Tj 



3. Optimal binning of the BABAR 2008 model 



Following Ref. [TH] the ratio of sensitivity to 7/^3 of a binned compared to an unbinned 
method is given by 



Q 2 
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(14) 



where 



Here, fs- is expressed as 



r i = J^fB-^dm^dm 2 ^ . 



Ib- = fD(m 2 + , m 2 J) + (x + iy)f D (m 2 _,m 2 + ) 



(15) 



(16) 



where x = r B cos (8b — 7) and y = tb sin (5b — 7). The parameter, Q, is the ratio of the 
number of standard deviations difference in the yields of B~ — > K~D° events as x and y 
change in a finite number of bins with respect to an infinite number of bins. The sensitivity 
to x and y is largely independent of their values. Therefore, again following Ref. [19], Eq. (14) 
can be simplified assuming x — y — to become 



Q \x=y=Q 



J2i Ni 



(17) 



where iVj is the number of B~ — > K~ D°(Kg7i + 7i^) events in the i th bin when is zero. 
Recalling that q and Sj are the amplitude-weighted averages of cosA<5d and sinA^ over 
each bin, respectively, it is clear that regions of similar A5e> will yield reasonable, though 
not necessarily optimal, values of Q. The Q value of the equal A5d binning presented in 



Sec. HI B 2 is 0.786 indicating that this binning choice is over 20% less sensitive statistically 
than an unbinned approach. (The values of Q are also computed for the different D° — > 
KgK + K~ binnings reported in Sec. |IIIAt the values are 0.771, 0.803, and 0.822 for Af = 2, 
Af 



3, and Af = 4 equal ASp. binnings, respectively.) 
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In Ref. [19] it was shown that the binning can be optimized to increase the value of Q. 
The optimization algorithm was provided by the authors of Ref. [19] and was adapted to 
use the lookup table of the BABAR 2008 model. The optimization is iterative and starts 



from the equal ASjj binning presented in Sec. IIIB2 Each iteration starts with the random 



selection of a sub-bin from the lookup table. In 90% of iterations the sub-bin is first tested 
to see if it lies on the boundary of a bin. If the sub-bin is not at a boundary the next 
iteration begins. Otherwise the sub-bin is moved from its current assignment to that of the 
neighboring bin and the value of Q is computed with the new assignment. If the value of Q 
is increased by this migration, the new assignment for this bin is kept and the next iteration 
begins. If the value of Q does not increase the assignment reverts to that originally given and 
the next iteration begins. In 10% of iterations, the selected sub-bin is given an assignment 
at random, irrespective of whether it is on a bin boundary. Again the reassignment of the 
sub-bin is kept if there is an improvement in Q; this allows part of one bin to 'grow' inside 
another bin if there is an improvement in the sensitivity. The procedure terminates when 
no further significant increase in the value of Q can be found. 

The binning that results from this optimization procedure is shown in Fig. |3](a) and is 
significantly different from that of the equal ASd binning (Fig. [T]). The optimized Q value is 
0.892, which is a 13% relative increase in sensitivity. However, there are many structures that 
are only a few sub-bins in size. Such regions are smaller than the experimental resolution 
and may result in systematic effects related to asymmetric migration of events from one bin 
to another. Furthermore, the position and shape of this fine structure depends critically 
on the components in the model, which may be realized differently in nature. Therefore, a 
smoothing procedure is implemented to remove these structures. The smoothing procedure 
starts by defining an 11 x 11 square of sub-bins centred about the sub-bin that is being tested. 
The number of sub-bins with the same bin assignment as the central sub-bin within the 
square is found. If the fraction of sub-bins of the same assignment is less than 30% the sub- 
bin assignment is changed to the modal bin assignment within the square. This procedure 
is performed for each sub-bin with the bin assignments from the original optimization to 
prevent bias. Figure [3](b) shows the binning after this smoothing procedure; the value of Q 
only decreases by 0.005. The smoothed optimal binning is used to calculate cp and sp . 

The same optimization procedure is applied to the three D° — > K®K + K~ binnings in 
terms of equal intervals of A<5q. However, the improvements in Q are found to be negligible 
compared to the equal binnings. Therefore, the cp and parameters for D° — > K^K + K~ 
decay are not reported here for optimized binnings. 



4- Modified- optimal binning of the BABAR 2008 model for the presence of background 

The Q values for the binnings provided are computed assuming that there is no back- 
ground present. There is a clear advantage to using the optimal binning in such a case, 
and simulation studies of a measurement of 7 using the observed values of q, Si, and the 
number of flavor-tagged D° — > K^-k + -k~ events in each bin, Ki, have confirmed the improved 
sensitivity in comparison to the equal A5d binning. However, when background is added 
to the simulation studies the sensitivity to 7 using the optimal binning can be worse than 



that for the equal A5n binning (see Sec. VII). The addition of background events naturally 



reduces the sensitivity to 7. The measurement of 7 is most sensitive when there are signifi- 
cant differences between yields in the bins for positive and negative B decays. In simulations 
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the 5-statistics according to the BABAR 2008 model, (b) The same binning after the smoothing 
procedure described in the text has been applied. The color scale represents the value 



there are two observed effects that can dilute the sensitivity. For the optimal binning there 
are bins where the asymmetry is large while the expected yields are low. If there is a large 
background yield in such a bin then the size of the asymmetry can be diluted to the point 
where the sensitivity gained by the optimal binning choice is significantly reduced by the 
presence of background. Also, very large bins contain a large fraction of the combinatoric 
background, which follows a reasonably uniform distribution over the Dalitz plot, which 
dilutes the asymmetry in that bin. With an assumed background model, it is possible to 
find a binning choice that maximizes the sensitivity to 7 in the presence of background. 

The background model assumed is determined from simulation studies of LHCb described 
in Ref. |16j . In this work three distinct types of background are considered. The first type of 
background is pure combinatoric, where the D° is reconstructed from a random combination 
of pions; the background-to-signal ratio, B/S, is expected to be less than 1.1 at the 90% 
confidence level. The second type of background is where a D meson is reconstructed 
correctly, and is then subsequently combined with a random kaon candidate to form a 
B candidate. The reconstruction of the D° and D is approximately equally likely and 
hence the distribution of this type of background in the i th bin will be proportional to 
(Ki + K_i). For this type of background, B/S is expected to be 0.35 ± 0.03. The third 
type of background, which has the smallest contribution to the total background, involves 
real B decays, predominantly B~ — > tt~ D° where the pion is misidentified as a kaon. In 
total the real B background has B/S less than 0.24 at the 90% confidence level. Sensitivity 
studies have shown that this type of background causes only a minor degradation in the 
sensitivity to 7. Therefore, this background type is ignored. In summary, the data sample 
is assumed to be composed of 41% signal events, 45% combinatoric background and 14% 
fully-reconstructed D background. 

In the presence of background the calculation of Q changes and will be written as Q' to 
distinguish it from the no background case. The value of Q' is still related to the number 
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of standard deviations by which the number of events in each bin is changed by varying 
parameters x and y, to the number of standard deviations if the Dalitz plot is divided into 



infinitely small regions as defined in Eq. (14), however the definition of is now 



\f B -? = f s \f D {m 2 + ,m 2 _) + {x + iy)f D (ml 



771 



\ 2 + f 1 B 1 (m 2 + ,m 2 _) + f 2 B 2 (mi. in: 



(18) 



where BUm 



+ ■ 



m 



and B 2 \n nr_ 



are the probability density functions for the combina- 



toric and fully-reconstructed D backgrounds, respectively, and f s , fx, and f 2 are the fractions 
of signal, combinatoric background, and fully-reconstructed D backgrounds, respectively. 
The assumed values of f s , fx and f 2 are 0.41, 0.45, and 0.14. 

As before the precision of x and y weakly depends on their values, therefore, the simplifi- 
cation that x = y = is once more made. In this case the expression analogous to Eq. (17) 
is given by 



E 



Q 



IJi + fxB Xi + f 2 B 



2 1 



--y=0 



PAf D {m 2 m 2 _)\ 2 \f D (m 2 _ } ml)\< 



f s \fD(ml,mt 



dm 2 , dm 2 



(19) 



I 2 + fxBx + f 2 B 2 

where Bu (B 2 i) is the integrated probability density functions for the combinatoric (fully- 
reconstructed D) background over the i th bin. 

The optimization algorithm to find the modified-optimal binning with the highest Q' is 
the same as described in Sec. Ill B 3 The modified-optimal binning Q' value is 0.910. In 
comparison, the equal ASu binning has Q' = 0.882 and the optimal binning has Q' = 0.867. 
The fine structure of the binning is smoothed out using the same technique as described 
for the optimal binning; the Q value drops by 0.006. The binning after the smoothing 
procedure is given in Fig. |4| In addition, we performed studies that show this binning choice 
retains the highest values of Q' even when the assumptions of the background model are 
modified. The alternative background models tested contain combinatoric background with 
a B/S between 0.8 and 1.3, and fully-reconstructed D° background with a B/S between 0.26 
and 0.44. 
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FIG. 5: Equal A5d binning of the D° — >K^-K + it~ Dalitz plot with M = 8 based on the Belle model 
|13j . The color scale represents the absolute value of the bin number, \i\. 

TABLE I: Comparison of the figure of merit values (Q, Q') calculated using the different models. 



Binning Fi£ 


Sure of merit BABAR 2008 p2] 


BABAR 2010 [14 


J Belle p3] 


Equal ASd 


Q 


0.786 


0.780 


0.762 


Optimal 


Q 


0.887 


0.880 


0.857 


Modified optimal 


Q' 


0.904 


0.903 


0.886 


Belle 


Q 


0.754 


0.755 


0.773 




5. Belle model binning 

The shape of the bins are dependent on the details of the amplitude model for the decay. 
There is another model available from the Belle experiment with which to define the bins 
[13]. This model does not use the same descriptions of the tttt and Kir S-wave as the BABAR 
2008 and 2010 models [T2"| IT4"]. A lookup table of this model has been provided by the Belle 
collaboration [2S]. Therefore, for a completeness, and to cross check our baseline results 
derived from the BABAR 2008 model, an equal A5d interval binning is derived from the 
latest Belle model. The binning over the Dalitz space is given in Fig. |5j 



6. Bin choice comparisons and the BABAR 2010 model 

The 2010 BABAR model only became available [24] after the completion of the D° — > 
Kg-K + 7T~ analysis. Therefore, to determine the possible impact on 7/^3 precision related 
to differences between the 2008 and 2010 models we compute the values of the respective 
figure of merit, Q or Q', for each BABAR model. This is done with the binnings fixed to 
those described above, which are derived from the BABAR 2008 model or the Belle model. 
The resulting figures of merit are given in Table [I] for the 2008 and 2010 BABAR models; 
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the values of are also given when computed with the Belle model. The values of 
computed using the BABAR 2010 model are slightly smaller than those for the 2008 model, 
but the difference is never greater than 0.007. In comparison, when calculating with 
the Belle model, for binnings derived from the BABAR 2008 model, the decrease is between 
0.018 and 0.030. Therefore, we conclude that using the BABAR 2008 model, rather than the 
BABAR 2010 model, to derive the binnings in this analysis will not result in a significant 
degradation in sensitivity to 7/03- 



IV. EVENT SELECTION 



This section summarizes the event selection for the two analyses. Section |IV A| describes 

K~ events. Section 



the selection of D° — > K SL K 



the selection of D 



IV B 



7T 



briefly summarizes the changes to 



events with respect to the previous analysis [7]. 



A. Selection of D° K% L K+K~ 

We perform the analysis on e + e~ collision data produced by the Cornell Electron Storage 
Ring at a center-of-mass energy, E cm , of 3.77 GeV. The data were collected by the CLEO-c 
detector and correspond to an integrated luminosity of 818 pb _1 . The CLEO-c detector 
is a solenoidal detector which includes a gaseous tracking system for the measurement of 
charged particle momenta and ionization energy loss, a ring-imaging Cherenkov detector 
to aid in particle identification, and a Csl crystal calorimeter to measure the energy of 
electromagnetic showers. The CLEO-c detector is described in detail elsewhere [30] . 

Samples of Monte Carlo (MC) simulated data are used to develop selection criteria, 
determine selection efficiencies, and to estimate certain types of background. EVTGEN [31] 
is used to generate the decays and GEANT |32j is used to simulate the CLEO-c detector 
response. Efficiency estimations are made on samples of signal events generated according 
to the D° — > KgK + K~ resonance model reported in Ref. [T2]. Separate signal samples are 
generated for each exclusive final state considered in the analysis and comprise 40,000 events 
per final state. Quantum correlations in the D°D° system are also simulated for each tag 
mode; this is particularly important for the CP-tagged D° decays. In addition a sample 
of generic DD decays corresponding to an integrated luminosity approximately 25 times 
greater than the data is used to estimate backgrounds. Quantum correlations are accounted 
for in the generic simulation. 

We adopt standard CLEO-c selection criteria for 7r + , tt°, and Kg mesons, which are 
described in Ref. [33]. The standard CLEO-c K + selection [33J is used for all final states 
apart from D° — > K s L K + K~ . For this final state, which has much smaller yields than 
K s L 7r + 7i~ , the significance of the signal is found to increase if the impact parameter criteria 
are loosened by a factor of four and the requirement on the fraction of associated tracking 
chamber hits compared to the expectation is removed. We require candidate K s — > tt + tt~ 
decays to have a mass within 7.5 MeV/c 2 of the nominal mass and the K s decay vertex is 
required to be separated from the interaction region by at least half a standard deviation. 
We reconstruct 77—7-77 candidates in a similar fashion to 7r° — > 77 candidates, with the 
requirement that the invariant mass is within 42 MeV/c 2 of the nominal mass; the same 
requirement is applied to r\ — )■ 7r + 7r~7r° candidates. Candidates for oj — > 7r + 7r~7r decays are 
required to be within 20 MeV/c 2 of the nominal u mass. We require rf — > 777r + 7r _ candidates 
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TABLE II: Single-tag (ST) and D° -> K° S L h + h~ double-tag (DT) yields. The single tag yields 
and uncertainties are computed following the method reported in Ref. [7] and are not corrected for 
efficiency The DT yields are the observed number of events in the signal region prior to background 
subtraction and before efficiency correction. 



Mode 


ST yield 




DT yields 








KgTT^TT 




K a s K + K~ 


K L K + K~ 


Flavor tags 












V— + 


"1 A A T* f O 1 A 

144563 ± 403 


1 A A A 

1444 


2857 


168 


302 


K 7V 7V 


258938 ± 581 


2759 


5133 


330 


585 


TS— 4- -+- — 


220831 ± 541 


2240 


4100 


248 


287 


K e v 




1 i m 

iiyi 




100 




Cr-even tags 












K + K~ 


13349 ± 128 


124 


357 


12 


32 


TV 7T 


6177 ± 114 


61 


184 


4 


13 




6838 ± 134 


56 




7 


14 






237 




17 










4 




T^O ^1 „+ „Q\ 
KjT]\T{ t 7T 7T J 








1 












4 












i 
i 




CP-odd tae-s 












K° S 7T 


19753 ± 153 


189 


288 


18 


43 




2886 ± 71 


39 


43 


4 


6 


K^t](tt + tt-tt°) 








2 


1 


K° s uj 


8830 ± 110 


83 




14 


10 


K° sV ' 

K° L 7T°7T 








3 
5 


4 






473 


1201 


56 


126 


K Q L -K + TT- 








140 




K° S K+K- 








4 
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to have an invariant mass in the range 950 to 964 MeV/c 2 . All nominal masses are taken 
from Ref. [23]. 

We consider Kg L K + K~ candidates reconstructed against the different final states listed 
in Table [IT} These are referred to as double-tagged (DT) events. More CP-tag final states are 
used in the analysis of K SL K + K~ than the K^ l tt + ii~ analysis [7] to increase the statistics 
available to determine q for this decay. (These modes are not included in the analysis of 
K^-k + -k~ because in this measurement the principal statistical limitation is the number of 
Kg-n + n^ vs. K1 s tt + tt~ events used to determine Sj.) We do not reconstruct final states 
containing two missing particles, such as K Q L K + K~ vs. K^tt . 

Final states that do not contain a K\ meson or neutrino are fully reconstructed via two 

kinematic variables: the beam-constrained candidate mass, = y / _E^ m /(4c 4 ) — p^/c 2 , 
where is the D candidate momentum, and AE = Ed — E cm /2, where Ed is the sum of 
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FIG. 6: (a) Distribution of mi, c for D° — > KgK + K candidates (D\) against the m& c for flavor-tag 
candidates (-D2): D° — >• K + ir~, K + ir~7r° and i^ + 7r _ 7r _ 7r + . The square signal region (red online) 



and four sideband regions are shown. Distributions of nib c for D° 
by (b) flavor, (c) CP-even, and (d) CP-odd decays. 



KqK + K candidates tagged 



the D daughter candidate energies. Signal decays will peak at the nominal D° mass and zero 
in rribc and AE, respectively. Mode-dependent requirements are placed on K° S K + K~ and 
tag candidates such that AE is less than three standard deviations from zero. The DT yield 
is determined from counting events in signal and sideband regions of the (mb c (D ), m^ c (I) )) 
plane in a manner similar to that presented in Refs. [2TI 131] . The signal region is defined 
as 1.86 GeV/c 2 < m bc (D°) < 1.87 GeV/c 2 and 1.86 GeV/c 2 < m bc (D°) < 1.87 GeV/c 2 . An 
example of the two-dimensional distribution of (mb c (D°), mi, c {D )) is shown in Fig. 6[a) for 
KgK + K~ candidates reconstructed against K + n~ , K + 7r~7r°, and K + 7t~7t~tt + decays. The 
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four different sidebands contain contributions from distinct types of combinatorial back- 
ground. The yields in these sidebands are scaled and subtracted from the yield in the signal 
region. The distributions for D° — > KgK + K~ candidates tagged by flavor, CP-even, 
and CP-odd final states are shown in Figs. |6](b), [6^c), and[6^d), respectively. These figures 
show clearly that the combinatorial backgrounds are small. The background-to-signal ratio 
for combinatoric background is less than 7.3% for all modes. 

To identify final states containing a single K\ meson, we compute the missing-mass 
squared recoiling against the fully-reconstructed D candidate and the particles from the 
other D decay containing the K\ meson. We select events consistent with the mass of the K\ 
meson squared. This technique was introduced in Ref. [35]. We reject events with additional 
charged tracks, 7r°, and rj candidates that are unassigned to the final state of interest. In 
addition, requirements are placed on any calorimeter energy deposits not associated with the 
charged or neutral particles that make up the final state of interest. The angle, a, between 
each unassigned shower and the missing-momentum direction is computed. Criteria are 
chosen to maximize signal sensitivity based on simulated samples of signal and background 
events. We retain events where cos a > 0.98, which indicates that the deposit is likely to be 
due to the interaction of the K\ meson with the calorimeter. When cos a < 0.98 mode-by- 
mode requirements are placed on the unassociated shower energy. The unassociated shower 
energy is required to be below a certain value which varies from 200 MeV for D° — > K^tt ^ 
candidates to 370 MeV for D° — > K\bj and D° — > K^rf candidates. Finally, criteria are 
placed on the momenta of 7r° and 77 candidates in tags containing a K\ meson to reduce 
background further. The combinatoric background yield in the signal region is estimated 
from the population in the lower and upper missing-mass squared sidebands. Information 
from the generic background simulation is used to determine the relative composition of the 
sidebands and the signal region to estimate better the combinatorial background. Figure [7] 
is the distribution of missing-mass squared for CP-tagged D° — > KgK + K~ candidates for 
data and simulated background, where the CP-eigenstate used to tag the event contains a 
K\ meson. 

We reconstruct the final state D° — > K + e~v by fully reconstructing a D° — > K®K + K~ 
candidate and requiring that the rest of the event contains both a kaon and an electron 
candidate of opposite charge. The quantity U miss = E miss — c|p miss | is used as a discrimi- 
nating variable, where P m i SS and p m i SS are the missing energy and momentum in the event, 
determined using the momenta of the fully reconstructed particles. The neutrino is the only 
particle not detected, so for a correctly reconstructed event U m i SS will equal zero. Figure [8] 
is the distribution of f/ m ; ss in the data and simulated background. The DT event yields for 
all final states are given in Table [TT| 

There are backgrounds to the signal that peak in tubc an d missing-mass squared at 
the same values as the signal, which can not be evaluated by examining the sidebands. 
These peaking backgrounds are estimated from the generic DD MC data samples. The 
largest peaking background to D° — > KgK + K~ decays is from D° — > K + K~h + 'k~ decays 
where the tt + tt~ pair form a Kg candidate. The peaking background from this source is 
estimated to be approximately 3.2% of the signal. For D° — > K < [K + K^ decays there are two 
significant peaking background contributions. The first source is D° — > K S K + K~ decays 
where the K$ is not reconstructed, usually in the case where it decays to 7r°7r°. The second 
significant source is from D° — > K + K^7i -k° decays where both the 7r° daughters are not 
reconstructed and the missing mass corresponds to that of the K\ meson. This background 
is not strictly peaking in that the 7r°7r° invariant mass does not always correspond to that 
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FIG. 7: Missing-mass squared distribution for CP-tagged D° — > K S K + K~ candidates where the 
CP-eigenstate contains a meson. The points are data and the solid histogram is the background 
estimated from simulation. 



of the K\ meson. However, the missing- mass squared distribution for D° — > K + K~ 7r°7r° 
background events determined from simulated data is not distributed linearly in the signal 
region and low missing-mass squared sideband, which would lead to a biased estimate of 
the background level if the sideband is used to determine the background level in the signal 
region. Therefore, the absolute level of this background is determined from the simulation 
and subtracted from both the signal and low missing-mass squared sideband. The D° — > 
K S K + K~ and D° — > K + K^tt tt peaking backgrounds are estimated to be 6.7% and 4.4% 
of the D° — > K^K + K~ signal, respectively. Tag modes that contain a K®, D — > K®X, 
also have a significant peaking background from D° — > KgX decays where the K s is not 



o 



70. 



reconstructed. The peaking background to D — > K\X events is between 4.0% and 6. 
The estimated peaking backgrounds are subtracted from the measured yields. The differing 
CP eigenvalues of the D K° S X and D ->• K° L X tags means that the distributions 
over Dalitz space of the signal and background can be significantly different. The effect of 
differing distributions of background is treated as a systematic uncertainty; the procedure 
to evaluate the uncertainty is described in Sec. VI B The background to K + e~v events is 
estimated as 1.8% in the signal region, denned as U m i SS < 50 MeV. 

We apply a kinematic fit to determine more reliably the position of a candidate in the 
Dalitz space. For final states containing a D° — > K$K + K~ decay, the fit constrains the 
invariant mass of both the signal and tag D° meson candidates to be the nominal D° mass 
and the Kg daughters to the nominal Kg mass. For D° — > K° L K + K~ decays there are 
several stages to the fit. The first stage constrains both the D° daughter kaons in the 
D° — > K < [K + K~ decay to originate from a common vertex and the tag decay candidate to 
the D° mass. The second stage uses the resulting four-momenta to estimate the mass of the 
missing K\ meson. The energy of the K\ candidate is rescaled such that the invariant mass 
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FIG. 8: U m i ss distribution for K + e~ ^-tagged K^K + K~ candidates. The points are data and 
the solid histogram is the background estimated from simulation. The background histogram (red 
online) is barely visible indicating that the background is negligible. 



is the nominal K\ meson mass. In the final stage the K\ candidate and the daughter kaon 
pair are constrained to the nominal D° meson mass. The introduction of the kinematic fit 
improves the resolution on m 2 + and m 2 _ by up to a factor of three for fully reconstructed 
KgK + K~ DT candidates. Although the fit only gives small improvements in core resolution 
for K\ final states there is a significant reduction in the number of events that lie in the 
non-Gaussian tails of the resolution distribution. 

The kinematic fit fails to converge in 1% to 3% of events, depending on the decay mode. 
In these cases the measured values of m 2 + and m 2 _ are rescaled such that they give the 
nominal D° mass. Even after the fit a small fraction of events are reconstructed outside the 
kinematic boundaries of the Dalitz space. The amplitude model used to assign events to a 
bin is undefined outside the Dalitz space. Therefore, the invariant masses of the candidate 
are changed to correspond to the point in the Dalitz space that is closest to the measured 
value in terms of the sum of (Am^_) 2 and (Am?.) 2 , where Am 2 is the residual between the 
measured value and a point within the Dalitz plot. 

The distributions of events across the Dalitz plane and as a function of m 2 K+K - for 
KgK + K~ reconstructed against CP-even tags are shown in Fig. [9j In addition, Fig. [9 
shows the distributions across the Dalitz plane and as a function of Tnr K+K - for K^K + K~ 
candidates reconstructed against CP-odd tags. The equivalent distributions for K®K + K~ 



candidates reconstructed against a CP eigenstate are shown in Fig. [101 The m 2 K+K - dis- 
tribution of K° S K + K~ (KlK + K~) candidates tagged with CP -even (CP-odd) eigenstates 
exhibits a peak due to the <fi resonance; as expected from CP conservation, this peak is 
not present for K§K + K~ (K®K + K~) candidates tagged with CP-odd (CP-even) eigen- 
states. Figure 11 shows the distribution of events across the Dalitz plane and as a function 
of 

m2 K+K- f° r K° S K + K candidates reconstructed against Kgn + ir decays. Furthermore, 
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FIG. 9: (a) Dalitz plot and (b) m 2 K+K _ distributions of D° — > K^K + K~ candidates tagged by a 
CP-even eigenstate. (c) Dalitz plot and (d) w? K+K _ distributions of D° — > KgK + K~ candidates 
tagged by a CP-odd eigenstate. 



shows the distribution of the same events across the K$tt 



Fig. [IT 

a function of the 7r + 7r~ invariant-mass squared, 
from using events in which both D mesons decay to Kg L h + h 
of flavor-tagged D° — > 



K° S K+K- 



and D 



and as a function of the ttik+k- are shown in Fig. [12 
to determine for each Dalitz-plot bin. 



7T _ Dalitz plane and as 
The increased statistics available 
is clear. The distributions 
K~ candidates across the Dalitz plane 
The flavor-tagged samples are used 
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FIG. 10: (a) Dalitz plot and (b) m 2 K+K _ distributions of D° — > K^K + K~ candidates tagged by a 
CP-even eigenstate. (c) Dalitz plot and (d) m 2 K+K - distributions of D° — > K^K + K~ candidates 
tagged by a CP-odd eigenstate. 



The efficiency in each bin is evaluated from the signal MC for each individual tag mode. 
Table III gives the total efficiency for each tag mode; these vary between (0.51 ± 0.04)% for 
K° S K + K~ vs. K° L n 7r° to (29.4 ± 0.3)% for K° L K + K~ vs. tt+tt". The uncertainty on the 
efficiency is that due to MC statistics. 

The finite detector resolution causes events to migrate between Dalitz-plane bins after 
reconstruction. Occasionally there is a significant asymmetric migration from one bin to 
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0.5 1.0 1.5 2.0 2.5 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2.0 2.2 



171+ (GeV 2 /c 4 ) m VtT) (GeV 2 /c 4 ) 

FIG. 11: (a) Dalitz plot and (b) m 2 K+K _ distributions of D° — > K®K + K~ candidates tagged 
by D -> K° s tt + tt- decays, (c) Dalitz plot and (d) mt + distributions of the D° — > K^tt + -k 
candidates in the same events. 

another. The effect is more pronounced in K^ L K + K~ decays than in Kg L 7T + n~ decays 
because of the presence of a relatively narrow and densely populated bin that encloses the 
<fi — > K + K~ resonance. We correct for this migration using MC data to determine the size 
and nature of the effect. For each binning we define a 2J\f x 2j\f matrix U for each DT mode 
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FIG. 12: (a) Dalitz plot and (b) m 2 K+K - distributions of flavor-tagged D — > KgK + K candidates, 
(c) Dalitz plot and (d) m 2 K+K _ distributions of flavor-tagged D° — > K^K + K~ candidates. 



as follows 

U» - , (20) 

k=-Af, k^0 

where rrij i is the number of signal MC events that are generated in bin j and reconstructed 
in bin i. The vector of migration-corrected yields in each bin, D corr , is determined from the 
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TABLE III: Selection efficiency for the different DT K SL K + K modes. The uncertainty is that 
due to MC simulation statistics. 



Tag 


Fffiripnrv (°7r\ 

n/iiicieiicy \/o) 




— 


K+ K~ 


TC~ 


r lavor tags 






l\ 7T 


i/i r -A- n 9 

14. D ± U.Z 


9^ 9 4- H q 
ZD.Z ± U.O 


iv 71 7T 


O.O It u.z 


1/1 q _|_ n 9 




in s + n 9 
lu.o in u.z 


k qxn 9 
u.y zn u.z 








CP-even tags 






K + K~ 


12.2 ± 0.2 


23 7 ± 3 




15.1 ± 0.2 


29.4 ± 0.3 




2.8 ± 0.1 


5.9 ± 0.1 




8.0 ± 0.1 




if? 77(77) 


7.9 ± 0.1 




1.6 ± 0.1 




Aj 

if? 7/ 


3.1 ± 0.1 




1.7 ±0.1 




G-r-OCLCL tagS 






s 


7 1-1-01 
/ . 1 zt U. 1 


in (i + n 9 

1U.U It u.z 




6.5 ±0.1 


9.7 ±0.2 


K° s 7](ir + TT-TT ) 


4.4 ±0.1 


6.8 ±0.1 


K° s uj 


3.4 ±0.1 


5.0 ±0.1 


K° S V' 


1.4 ±0.1 


2.0 ±0.1 


K L 7T 7T 


0.51 ±0.04 




K° s ir + ir- 


7.9 ±0.1 


13.0 ±0.2 




14.0 ±0.2 




K° S K+K- 


4.9 ±0.1 


7.0 ±0.1 



vector of reconstructed yields in each bin, D rec , using the relation: 

D corr = U-l D rec_ ^ 

As an example the migration matrix for K Q L K + K~ vs. K^tt events when the D° — y 
K$K + K~ Dalitz space is divided into M = 3 bins is given in Table IV Typically the 



migration out of the bin containing the <fi resonance is 5% for D° — y K$K + K~ modes and 
between 10% to 25% for D° — y K^K + K~ modes. The errors on the elements of U due to 
the limited MC statistics are treated as a systematic uncertainty. 



B. Updates to the Kg L 7T + 7r selection 

Single tags are used in the analysis of D° — y K^-k + -k^ to determine the normalization 



factors in Eqs. fej) and (12) with limited systematic uncertainties. (In the D — y KgK + K 



analysis these systematic uncertainties are not as important given the available statistics and 
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TABLE IV: Migration matrix U elements (%) for K^K + K vs. K s tt° events when the D° — y 
K^K + K~ Dalitz space is divided into M = 3 bins. 



?' 


U it i 


U it 2 


Ui, 3 


Ui,-i 


U it - 2 


Ui,- 3 


1 


86.2 


11.8 


0.5 


0.4 


0.0 


0.0 


2 


11.3 


88.1 


0.5 


0.0 


0.0 


0.0 


3 


1.5 


0.1 


98.9 


0.0 


0.0 


0.0 


-1 


1.0 


0.1 


0.0 


85.2 


10.6 


0.0 


-2 


0.0 


0.0 


0.0 


13.9 


89.3 


1.0 


-3 


0.0 


0.0 


0.0 


0.4 


0.1 


99.0 



therefore the normalizations are determined from the number of DD pairs and the measured 
branching fractions.) The ST selection is identical to that described in Ref. [7J. However, 
an error in the previous analysis led to a small fraction (3.6%) of the data being excluded 
from the ST analysis of CP eigenstates. The updated yields are given in Table |H] 

The selection requirements of DTs are identical to that in the previous analysis [7J. 
However, there are several changes to the yields reported. The same data excluded from 
the original ST analysis were missing from the evaluation of yields of final states containing 
a K\ candidate. 2 Also, the kinematic fit procedure described in Sec. IV A has also been 
applied to the i^°7r + 7r~ events. However, events that lie outside the Dalitz plot after the fit 
procedure are rejected rather than migrated into the physical region, which changes some 
of the yields compared to those reported in Ref. [TJ; in the previous analysis all events were 
retained. The updated DT yields are shown in Table [TT| 



V. 



EXTRACTION OF a AND s< 



The efficiency-corrected, background-subtracted, and migration-corrected bin yields for 



each DT mode need to be normalized to determine the measured values of Af/'^ and so 
that they can be used to evaluate c(' and sj via Eqs. (8 



r(') 



(10), (11), and (12). In addition, 



the values of flavor-tagged yields in each bin must also be normalized appropriately to obtain 
K\ ; . For the KqK + K~ analysis the single-tag yields in the normalization factors hep and 



s0 



in Eqs. (8), (10), (11), and (12) are determined from the number of D D pairs in 



the sample, N DQ jp , multiplied by the branching fractions of the modes taken from Ref. 
The value of N DQ -^o is calculated from the integrated luminosity and the value of the cross 

section for e + e~ — y ^(3770) — y D°D° reported in Ref. [33]. For the l ii + -k~ analysis the 
measured single-tag yields are used for all modes apart from those tagged by D° — > K^n and 
D° K~e + v, which cannot be reconstructed exclusively; in these cases the normalization 
is performed as in the D° — y Kg L K + K~ analysis. 

Yields of K s L h + h~ events selected against Cabibbo-favored (CF) hadronic flavor tags are 
contaminated with doubly-Cabibbo suppressed (DCS) decays. We refer to these hadronic 



2 In addition, the yields of K^ti + h 
in Table II of Ref . [7]. 



K e + v and K ( ii: + 'K vs. K%tt + t: were incorrectly documented 
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TABLE V: Values of the parameters used to make the corrections to the pseudo-flavor tag yields 
and the references from which they are taken. 

Ktt 5.80 ± 0.08 [37] 202 ± 10 [37] 1 

Kirir 4.8 ±0.2 [38] 227 ±17 [21] 0.84 ± 0.07 [2T] 

i^vrvrvr 5.7 ± 0.2 [38] 114 ± 26 |2T] 0.33 ± 0.26 [ZQ 



TABLE VI: Values of F/^ (%) measured from the flavor-tagged D° ->■ AT^AT+A' - data for the 
M = 2, M = 3, and A/" = 4 equal AJp binnings. Predicted values from the BABAR 2010 model of 
D° — > KgK + K~ are also given. 



i 


F (%) 


(%) 




Measured 


Predicted Measured 


Predicted 




N 


= 2 equal A#d binning 




1 


23.9 ± 1.6 


22.5 ±4.2 35.5 ±1.9 


28.6 ±1.1 


2 


17.3 ± 1.5 


21.1 ± 1.2 23.3 ±1.7 


27.8 ±4.1 




N 


= 3 equal A8d binning 




1 


22.0 ± 1.5 


19.8 ±3.8 33.0 ±1.7 


25.6 ± 1.0 


2 


18.1 ± 1.4 


22.7 ± 1.4 22.8 ±1.6 


26.1 ±5.3 


3 


1.2 ±0.4 


1.4 ±0.7 3.0 ±0.6 


3.8 ± 1.6 




A" 


= 4 equal A5o binning 




1 


20.0 ± 1.5 


18.3 ±3.3 30.5 ±1.7 


23.0 ± 1.1 


2 


7.2 ± 1.1 


8.5 ±1.0 7.6 ±1.3 


8.6 ± 1.3 


3 


13.3 ± 1.4 


16.3 ±1.3 17.7 ±1.4 


21.3 ±4.0 


4 


0.8 ±0.4 


0.5 ±0.4 2.8 ±0.6 


3.5 ± 1.3 



flavor modes as pseudo-flavor tags. This introduces a bias in the extracted values of 
[7J. To account for this effect, the flavor-tagged yields in each bin are scaled by a correction 
factor, which is estimated using the D° — > Kg L K + K~ and D° — > Kg l tt + tt~ decay models 
reported in Refs. [H] and [TT], respectively. The correction factor for pseudo- flavor tag, 
F G (K-7T + , K-7T+7T , ir-vr+vr+vr-), is 

/, \f(rn 2 i _,m 2 _)\ 2 dm 2 i _dm 2 _ 

Si{\f{rn? + ,ml)\ 2 + {r^) 2 \f{m 2 _, m 2 + )\ 2 + 2r^R F ^[e~ iS D f{m 2 + , m 2 _)f*{m 2 _ , m 2 + )])dm\dm 2 _ ' 

(22) 

where is the ratio of the DCS to CF decay amplitudes and 8^ is the associated average 
strong-phase difference. Rf is the coherence factor for decays to three or more particles [36] 
and equals unity for two-body decays. The values of these parameters and the references 
from which they are taken are given in Table |Vj 

Equation ([3| defines the normalized values of the flavor-tagged yields in each bin, 

such that = Kr^i/Ao, where A D = YjiLi{Ki + K_i). The fully-corrected values of 

Ft^i measured for the M = 2, Af = 3, and Af = 4 equal ASd binnings are given in Table VI 
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TABLE VII: Values of •?(_)» (%) measured from the flavor-tagged D° — > Kgir + ir data for the 
optimal binning. Predicted values from the BABAR 2008 model of D° — > K^tt + tt^ are also given. 



Ft (%) (%) 
Measured Predicted Measured Predicted 



1 


17.0 ±0.5 


17.2 


8.3 ±0.4 


8.3 


2 


8.4 ±0.4 


8.4 


2.4 ±0.2 


1.9 


3 


7.2 ±0.3 


6.9 


2.3 ±0.2 


2.0 


4 


2.4 ±0.2 


2.5 


1.6 ±0.2 


1.7 


5 


7.6 ±0.4 


8.6 


4.8 ±0.3 


5.3 


6 


5.9 ±0.3 


5.9 


1.3 ±0.2 


1.5 


7 


12.8 ±0.5 


12.4 


1.6 ±0.2 


1.4 


8 


13.0 ±0.5 


13.0 


3.1 ±0.2 


2.8 



TABLE VIII: Values of Fr ^ (%) measured from the flavor-tagged D° — > K^tt + tt~ data for the 
equal A5d binning derived from the BABAR 2008 model. Predicted values from the BABAR 2008 
model of D° — > Kgir + ir~ are also given. 

i Fi (%) F^ (%) 

Measured Predicted Measured Predicted 



1 


9.0 ±0.4 


9.1 


2.6 ±0.2 


2.2 


2 


14.4 ±0.5 


14.2 


0.5 ±0.1 


0.5 


3 


14.7 ±0.5 


14.1 


0.3 ±0.1 


0.4 


4 


9.9 ±0.4 


10.0 


5.9 ±0.3 


6.6 


5 


5.7 ±0.3 


5.4 


3.3 ±0.2 


3.2 


6 


7.5 ±0.4 


7.4 


0.5 ±0.1 


0.4 


7 


10.9 ±0.4 


11.5 


5.5 ±0.3 


5.4 


8 


2.2 ±0.2 


2.2 


6.9 ±0.3 


6.3 



The results are the average of the pseudo-flavor tag modes with the D° — > K~e + v tagged 
data. Also given are the predicted values for from the BABAR 2010 amplitude model 
of D° — > K$K + K~ decays. The error on the predicted value of Fim is determined from the 
uncertainties on the amplitude-model parameters. The agreement between the measured 
and predicted values is reasonable in all bins except F_i, which is different by more than 
three standard deviations for all binnings. This discrepancy may be a statistical fluctuation 
or indicate a deficiency in the model in this region. In order to ascertain whether this 
effect can lead to a significant bias in our analysis, we perform the fit to determine cf 1 and 
s( with the predicted rather than measured values of K^y, the difference in fit results is 
negligible. Therefore, we conclude it is reasonable to use our measured values of to 
determine the parameters. Tables VII, VIII, IX, and |X| show the measured values of for 
the D° — > Kgir + 7T~ data divided according to the optimal, BABAR 2008 model with equal 
A5d, modified-optimal, and Belle equal A5d binnings, respectively. Again these results are 
the average of pseudo-flavor and semileptonic tagged data. Predicted values are also given. 
In this case the uncertainty due to the amplitude-model parameters is negligible compared 
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TABLE IX: Values of i*(_v (%) measured from the flavor-tagged D° — > K®it +, k~ data for the 
modified-optimal binning. Predicted values from the BABAR 2008 model of D° — > Kgir^ir" are 
also given. 

~i Fi (%) (%) 

Measured Predicted Measured Predicted 



1 


5.4 ±0.3 


5.1 


1.5 ±0.2 


1.6 


2 


16.0 ±0.5 


16.2 


2.1 ±0.3 


2.0 


3 


22.0 ±0.6 


21.6 


2.3 ±0.3 


2.1 


4 


7.8 ±0.4 


8.8 


4.9 ±0.5 


5.3 


5 


3.8 ±0.3 


3.9 


3.1 ±0.4 


2.8 


6 


8.3 ±0.4 


8.1 


1.1 ±0.2 


1.2 


7 


8.7 ±0.4 


9.0 


4.4 ±0.5 


4.6 


8 


2.2 ±0.2 


2.3 


6.0 ±0.6 


5.5 



TABLE X: Values of (%) measured from the flavor-tagged D° — > KgTr + iT~ data for the equal 
A5d binning derived from the Belle model. Predicted values from the BABAR 2008 model of 
D° — > Kgir + ir~ are also given. 

~i Fi (%) F_i (%) 



Measured Predicted Measured Predicted 



1 


16.5 ±0.5 


16.5 


8.8 ±0.4 


8.0 


2 


7.7 ±0.4 


7.6 


2.0 ±0.2 


1.6 


3 


9.8 ±0.4 


10.2 


3.2 ±0.2 


2.8 


4 


3.0 ±0.2 


3.0 


1.3 ±0.1 


1.2 


5 


8.0 ±0.4 


9.2 


4.0 ±0.3 


4.6 


6 


7.1 ±0.3 


7.3 


1.8 ±0.2 


1.7 


7 


9.9 ±0.4 


10.0 


1.6 ±0.2 


1.3 


8 


12.4 ±0.4 


12.2 


2.9 ±0.2 


2.6 



to the uncertainties on the measurements. The agreement between measured and predicted 
values is reasonable in all cases. 

We use the corrected and normalized values of the bin yields to determine cp and s\ . 
The fits to the Kg L 7i + 7T^ and Kg L K + K~ data are made separately. We perform the 
fit to Kg L Tr + Tr~ data first because the Kg L K + K~ fit depends upon the values of the 
cp and s( for D° — > Kg L ir + TT~ decays when including K$K + K~ vs. K® L 7r + 7r~ and 
K° L K + K~ vs. Kg7c + iT~ candidates. The fit results from the equal A5d binning derived 
from the BABAR 2008 model are used in the fit to K® L K + K~ data; the Kg L 7i + n~ strong- 
phase difference parameters are fixed to the measured values in the nominal fit, these are 
then varied within their errors to determine the related systematic uncertainty (Sec. VI B). 
In the Kgir + TT~ analysis we obtain values of cp and sj by minimizing the log-likelihood 
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expression 



21og£ = -2^1ogP(M±,(M±)) K o 7r+7r - iCP 

i 

-2$>gP(M*, (Mf )) Kl ^-, CP 

i 

-2j2^gP(M ij: (My ) ) K°TT+TT-,K°-K+Tr 

-2X)logP(JWj,(Mi>W w 



X"°7r+7r- 



In the K$K + K analysis we minimize the expression 



(23) 



- 21og£ 



-2^1ogP(M±, (M±)) K o K+K - >CP 

i 

-2£logP(Mf , {M^)) nK+K - t cP 

i 

— 2^ logP(Mjj, (Mij}) K o K+K - K o K+h 
-2^]ogP{Mi d ,(M' 1 j)) K o K+K 
-2^ogP(M lv (M tJ )) K o K+K 

id 

-2^1ogP(M4,(M;.)) x o 

K+K 

id 

-2^ogP(M: p (M' tJ }) K o K+K 



K° L K+K- 



K Q S TI+TV- 



+x 2 - 



(24) 



Here the expected number of CP-tagged K®h + h~ (Klh + h~) events in the i th bin, (Mf) 
((M^)), is determined from Eq. (|8j) (Eq. (11)). Similarly the expected number of events 
where both D mesons decay to Kg L h + h~, (Mij) ((ML)) is determined using Eq. (10) 
(Eq. (12)). The function P(M, (M)) is the Poisson probability of obtaining M events given 
a mean of (M). There is an additional \ 2 term 



X 



E 



Ac, 



8Ac, 



+ E 



s ■ 



As, 



SAsi 



(25) 



which constrains the extracted c, (s'J to differ within errors from q (sj) by the predicted 
quantities Aq (As*). 

We briefly discuss the estimation of Ac, (As,) and their uncertainties. An isobar res- 
onance model must be used to determine these constraints; we use the D° — > K^K + K~ 
and D° — > KgTT + 7c~ models reported in Refs. [H] and [Uj, respectively. The intermediate 
resonances contributing to D° — > K^h + h~ model differ in two ways from those contributing 
to D° — > Kgh + h~ . Firstly, DCS decays contribute with opposite sign. This can be seen by 
considering the D° — > Kg L h + h~ amplitudes, A, in terms of those to the flavor eigenstates 
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D° -»■ K°h + hr and D° -»■ FV/r 

K° s h + h-) = ^= \A(D° ->■ K°h + h-) + ^(D ZV/T)1 (26) 

„4(D° -> K° L h + h-) = 1 \A(D° -> K°h+h-) - A(D° -> 3rV/T)l . (27) 

V2 

The relative minus sign between the terms in ^4(D° — >■ K^h + h~) can be accommodated 
by introducing a 180° phase difference for all DCS contributions to the _D° — > K®h + h~ 
model compared to the same contribution to the _D° — > Kgh + h~ model. Secondly, for CP- 
eigenstate amplitudes, such as D° — > Kg L <f>, the D° — > K^h+h^ amplitude can be related 
to the D° Kgh + h~ amplitude by multiplying the latter by a factor (1 — 2re 4<5 ), where 
r is of the order tan 2 8c and 5 can take any value. Here, 9c is the Cabibbo angle. The 
origin of this factor is again related to the sign difference between DCS contributions to 
the D° — > Kgh + h~ and D° — > K^h + h~ amplitudes, and is analogous to the mechanism 
which induces the difference in rates for D — > K^tt and D — > K\ii decay [35j [39]. We 
determine central values of Aq and As« by assuming r = tan 2 $c |23j and 5 = 0°. Part 
of the uncertainty on Aq and Asi is evaluated by randomly choosing the assumed values 
of r and 5 100 times, and recomputing the constraints for each set of parameters. The 
value of 5 is assumed to have a equal probability to lie between 0° and 360° and that of 
r to have a Gaussian distribution of mean tan 2 9 C and width 0.5 x tan 2 C . The RMS of 
the resulting distributions of Aq and Asj are taken as the uncertainties from this source. 
A second source of uncertainty is related to the model choice. For the D° — > K®K + K~ 
model this is estimated by varying the isobar model parameters by their uncertainties [H] , 
accounting for any correlations among the parameters [21], then recomputing Aq and As,. 
The differences with respect to the values of the constraints computed with the nominal 
values of the parameters are taken as the uncertainty on Aq and Asj from this source. For 
D° — > Kg7T + TT~ we follow Ref. [7] and consider two alternative models JTUJ HO] to determine 
Aq and Asf, the largest deviation of the central value from that computed with the default 
model [TT] is taken as a systematic uncertainty. As examples of the constraints found by this 
procedure, the values of Acj and Asj for the M = 3 division of the D° — > KgK + K~ Dalitz 
plot and the optimal binning of the D° — > f^7r + 7r~ Dalitz plot are given in Tables XI and 



XII , respectively. The size of Aq and As« can be significant for bins dominated by either 
a DCS decay, such as the i^* + (892) in bin three of the optimal D° — > Kg7t + 1T~ binning 
(Fig. |3^b)), or a neutral resonance, such as the a°(1450) in bin three of the J\f = 3 division 
of the D° -> K° S K + K~ Dalitz plot (Fig. Q. 

TABLE XI: Values of Aq and Asj constraints for the M = 3 equal ASd binning of the D° — > 
K° S K+K- Dalitz plot. 



i 


Aa 


As., 


1 


0.026 ±0.014 


-0.007 ±0.023 


2 


0.041 ±0.019 


0.012 ±0.010 


3 


-0.563 ±0.311 


0.713 ±0.161 



The contribution of the \ 2 t° the likelihood is investigated to ensure that these constraints 
are not leading to any significant bias. For the fits to D° — > Kg L K + K~ and D° — > K s L 7r + 7r~ 
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TABLE XII: Values of and Asj constraints for the optimal D — y i^^7r + 7r binning. 



i 


Aa 


Asi 


1 


0.39 ±0.17 


0.07 ±0.06 


Z 


U.lo It u.uo 


n m -i- n i n 


3 


0.61 ±0.15 


0.30 ±0.12 


4 


0.09 ±0.08 


0.00 ±0.08 


5 


0.16 ±0.17 


0.06 ±0.06 


6 


0.57 ±0.21 


-0.15 ±0.24 


7 


0.03 ±0.01 


-0.04 ±0.06 


8 


-0.10 ±0.15 


-0.15 ±0.21 



the ranges of x 2 are 0.30 to 0.75 and 1.0 to 2.3, respectively. In addition, no individual bin 
contributes more than 0.9 to the total x 2 ■ Therefore, we conclude that the constraint is 
not biasing our result significantly from the values favored by the data, and is principally 
improving the precision of the parameters. 

fitting procedure has been tested using samples of signal MC 



The D° -»■ 



events. In the validation procedure, the number of events for each tag is assumed to follow 
a Poisson distribution about the mean expectation, while the ratio between double-tagged 
events and single-tagged events is computed as B K o K+K -/2, ignoring quantum-correlations. 
The means of the fitted q and Sj distributions exhibit small, but statistically significant, 
biases due to the assumptions made in the fit. The systematic uncertainty we associate to 
the bias is described in Sec. IVI Bl 

The parameters that result from fits to the D° — y Kg L K + K~ data divided into Af = 2, 
M = 3, and M = 4 equal A5d bins are given in Tables XIII and XIV The statistical 
uncertainty on the parameters includes that related to the Aq and Asj constraints. The 
statistical correlations among the parameters for the M = 3 equal ASd binning are shown 
in Table XV. The other statistical correlations are given in Ref. [28J. 

give -> K° s 



Tables 



XVI 



and 



XVII 



data divided according 

to the four different binnlngs. The statistical correlation matrices for each binning are given 
in Ref. 



VI. SYSTEMATIC UNCERTAINTIES 



The systematic uncertainties on the measured values of cp and s\'' come from a variety of 



-,(') 



sources. Tables XVIII and XIX give examples of the individual sources of uncertainty for the 
M = 3 equal ASd binning of the D° — y K S K + K~ Dalitz plot and the equal A5d binning of 
the D° — y Kg7r + ir~ Dalitz plot, respectively. The breakdown of the systematic uncertainty 
is similar for the other binnings considered for D° —y K$K + K~ and D° — y K^h + 'k~ data. 
Many of the systematic uncertainties are from common sources for the two final-states; 



these are described in Sec. |VI A[ The sources of uncertainty considered exclusively for the 
D° — y KgK + K~ analysis, or those that are evaluated in a significantly different fashion than 
D° — y K%-k + -k^ , are described in Sec. |VI Bl A discussion of the uncertainties related to the 



'0 



7T 



analysis are given in Sec. VI C A brief 



acceptance and background to the D° — y K 
discussion of the uncertainty related to the constraints Aq and Asj is given in Sec. |VI D 
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TABLE XIII: Measured values of q and Si for the different D — > K$K + K~ binnings. The first 
uncertainty is statistical, including that related to the Aq and As; constraints, and the second 
uncertainty is systematic. 



% 


Ci 


Si 




Jx = I equal lAod 


bins 


1 


0.818 ± U.1U7 ± 0.037 


0.445 ± 0.215 ± 0.143 


o 
z 


— U. (46 ± 0.083 ± 0.035 


U.zzy ± U.zzU ± 0.0/y 




J\f = 3 equal A<5d 


bins 


1 


0.793 ± 0.063 ± 0.038 


0.431 ±0.222 ±0.142 


2 


-0.566 ± 0.092 ± 0.034 


0.413 ±0.234 ±0.094 


3 


-0.096 ±0.329 ±0.131 


-0.461 ± 0.432 ± 0.175 




Af = 4 equal A5 D 


bins 


1 


0.858 ± 0.059 ± 0.034 


0.309 ±0.248 ±0.180 


2 


0.176 ±0.223 ±0.091 


0.992 ±0.473 ±0.403 


3 


-0.819 ± 0.095 ± 0.045 


0.307 ±0.267 ±0.201 


4 


0.376 ±0.329 ±0.157 


-0.133 ±0.659 ±0.323 



TABLE XIV: Measured values of c- and s- for the different D° ->• K° S K + K~ binnings. The first 
uncertainty is statistical, including that related to the Aq and Asj constraints, and the second 
uncertainty is systematic. 



i 




s'- 




M = 2 equal AJp 


bins 


1 


0.839 ±0.108 ±0.073 


0.445 ±0.216 ±0.150 


2 


-0.775 ± 0.085 ± 0.068 


0.298 ±0.220 ±0.093 




M = 3 equal AJp 


bins 


1 


0.814 ±0.063 ±0.064 


0.422 ± 0.222 ±0.143 


2 


-0.529 ±0.092 ±0.071 


0.426 ±0.234 ±0.098 


3 


-0.583 ±0.349 ±0.197 


0.241 ±0.456 ±0.181 




= 4 equal ASd 


bins 


1 


0.874 ±0.059 ±0.055 


0.303 ±0.248 ±0.180 


2 


0.270 ± 0.225 ±0.160 


0.965 ±0.473 ±0.406 


3 


-0.810 ± 0.096 ± 0.060 


0.346 ±0.268 ±0.204 


4 


-0.317 ± 0.408 ± 0.201 


0.770 ±0.696 ±0.350 



A. Common systematic uncertainties 

The statistical uncertainties on the measurements of using flavor and pseudo-flavor 
tagged events are not included in the fit to determine q, Sj, c'j, and s^. Even though 
there are significantly more flavor and pseudo-flavor tagged events than the CP-tagged or 
Kg L h + h~ vs. K^h + h~ events, there is a non-negligible uncertainty related to the limited 
statistics of these samples. Each input is varied separately by its uncertainty and the fits are 
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TABLE XV: Statistical correlations (%) among the parameters for N = 3 equal A8d binning of 
the D° -> KgK + K~ Dalitz plot. 





C2 


C3 


si 


S2 


S3 


4 


4 


4 


s'i 


4 


4 


Cl 


0.6 


-2.9 


-1.3 


-0.6 


-0.2 


97.6 


0.6 


-2.1 


-1.3 


-0.6 


-0.1 


C2 




0.8 


1.1 


3.2 


-0.0 


0.6 


98.0 


0.4 


1.1 


3.2 


-0.0 


C3 






-0.3 


0.2 


2.5 


-2.8 


0.8 


66.4 


-0.3 


0.2 


2.4 


Si 








-2.0 


6.1 


-1.3 


1.1 


-0.1 


99.4 


-2.0 


5.8 


S2 










-3.3 


-0.5 


3.2 


0.1 


-1.9 


99.9 


-3.0 


S3 












-0.2 


-0.0 


2.2 


6.1 


-3.3 


93.7 


4 














0.6 


-2.0 


-1.3 


-0.5 


-0.1 


4 
















0.4 


1.1 


3.2 


-0.0 


4 


















-0.1 


0.1 


2.2 






















-1.9 


5.8 


4 






















-3.0 



TABLE XVI: Measured values of q and Sj for the 
uncertainty is statistical, including that related to 
uncertainty is systematic. 

Optimal 

Ci Si 

-0.009 ± 0.088 ± 0.094 -0.438 ± 0.184 ± 0.045 
0.900 ± 0.106 ± 0.082 -0.490 ± 0.295 ± 0.261 
0.292 ± 0.168 ± 0.139 -1.243 ± 0.341 ± 0.123 
-0.890 ± 0.041 ± 0.044 -0.119 ± 0.141 ± 0.038 
-0.208 ± 0.085 ± 0.080 0.853 ± 0.123 ± 0.035 
0.258 ± 0.155 ± 0.108 0.984 ± 0.357 ± 0.165 
0.869 ± 0.034 ± 0.033 -0.041 ± 0.132 ± 0.034 
0.798 ± 0.070 ± 0.047 -0.107 ± 0.240 ± 0.080 



different D° — > KgTT + -rr binnings. The first 
the Aa and Asi constraints, and the second 

Equal A5 D BABAR 2008 

C-i Si 

0.655 ± 0.036 ± 0.042 -0.025 ± 0.098 ± 0.043 
0.511 ± 0.068 ± 0.063 0.141 ± 0.183 ± 0.066 
0.024 ± 0.140 ± 0.080 1.111 ± 0.131 ± 0.044 
0.569 ± 0.118 ± 0.098 0.328 ± 0.202 ± 0.072 
0.903 ± 0.045 ± 0.042 -0.181 ± 0.131 ± 0.026 
0.616 ± 0.103 ± 0.072 -0.520 ± 0.196 ± 0.059 
0.100 ± 0.106 ± 0.124 -1.129 ± 0.120 ± 0.096 
0.422 ± 0.069 ± 0.075 -0.350 ± 0.151 ± 0.045 



Modified optimal 



Equal A5 D Belle 



Ci 



Si 



-0.216 ±0.104 ±0.088 
0.827 ±0.060 ±0.053 
0.101 ±0.085 ±0.118 
-0.955 ±0.038 ±0.034 
-0.522 ±0.095 ±0.079 
0.291 ±0.102 ±0.075 
0.682 ±0.051 ±0.047 
0.724 ±0.071 ±0.044 



-0.399 ±0.204 ±0.049 
-0.031 ±0.172 ±0.084 
-0.558 ±0.161 ±0.070 
-0.204 ±0.137 ±0.055 
0.911 ±0.130 ±0.067 
1.030 ±0.196 ±0.065 
-0.037 ±0.146 ±0.029 
-0.180 ±0.194 ±0.050 



0.710 ±0.034 ±0.038 
0.481 ±0.080 ±0.070 
0.008 ±0.080 ±0.087 
-0.757 ±0.099 ±0.065 
-0.884 ±0.056 ±0.054 
-0.462 ±0.100 ±0.082 
0.106 ±0.105 ±0.100 
0.365 ± 0.071 ± 0.078 



-0.013 ±0.097 ±0.031 
-0.147 ±0.177 ±0.107 
0.938 ±0.120 ±0.047 
0.386 ±0.208 ±0.067 
-0.162 ±0.130 ±0.041 
-0.616 ±0.188 ±0.052 
-1.063 ±0.174 ±0.066 
-0.179 ±0.166 ±0.048 
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TABLE XVII: Measured values of d and d for the different D c 



K° S TT + TT- 



binnings. The first 



uncertainty is statistical, including that related to the Ac, and Asj constraints, and the second 
uncertainty is systematic. 



Optimal 



Equal A5 D BABAR 2008 



s' 



0.470 
1.073 
0.869 
-0.786 
-0.139 
0.654 
0.901 
0.817 



± 0.096 
±0.102 
±0.142 
± 0.047 
± 0.089 
±0.135 
± 0.034 
± 0.090 



±0.102 
±0.126 
±0.165 
±0.052 
±0.093 
±0.152 
± 0.047 
± 0.095 



-0.363 ± 
-0.501 ± 
-0.890 ± 
-0.137 ± 
0.921 ± 
0.832 ± 
-0.076 ± 
-0.157 ± 



0.185 ±0.191 
0.297 ±0.397 
0.329 ±0.352 
0.154 ±0.159 
0.126 ±0.132 
0.326 ±0.362 
0.132 ±0.137 
0.281 ±0.306 



0.768 ± 
0.679 ± 
0.278 ± 
-0.446 ± 
-0.824 ± 
-0.240 ± 
0.480 ± 
0.708 ± 



0.038 ±0.051 
0.067 ±0.090 
0.106 ±0.112 
0.116 ±0.128 
0.051 ±0.056 
0.116 ±0.123 
0.106 ±0.113 
0.066 ± 0.083 



-0.079 ±0.095 
0.080 ±0.183 
1.090 ±0.106 
0.455 ±0.219 
-0.194 ±0.136 
-0.557 ±0.201 
-0.975 ±0.104 
-0.285 ±0.150 



±0.108 
±0.196 
±0.119 
±0.235 
±0.140 
±0.211 
±0.141 
±0.158 



Modified optimal 



Equal A5 D Belle 



s' 



-0.049 
0.935 
0.614 
-0.876 
-0.357 
0.584 
0.789 
0.717 



±0.108 
± 0.056 
± 0.082 
± 0.047 
± 0.094 
± 0.094 
± 0.058 
± 0.080 



±0.125 
±0.072 
±0.089 
±0.051 
±0.100 
±0.103 
±0.070 
± 0.089 



-0.386 ± 
-0.028 ± 
-0.398 ± 
-0.249 ± 
0.980 ± 
0.963 ± 
-0.091 ± 
-0.219 ± 



0.218 ±0.226 
0.168 ±0.191 
0.160 ±0.175 
0.139 ±0.149 
0.131 ±0.151 
0.200 ±0.214 
0.145 ±0.149 
0.201 ±0.209 



0.817 ± 
0.668 ± 
0.197 ± 
-0.592 ± 
-0.720 ± 
-0.121 ± 
0.426 ± 
0.641 ± 



0.035 ± 0.048 
0.079 ±0.104 
0.082 ± 0.090 
0.117 ±0.125 
0.056 ±0.062 
0.108 ±0.117 
0.104 ±0.113 
0.071 ± 0.089 



-0.070 ±0.095 
-0.219 ±0.176 
0.935 ±0.117 
0.520 ±0.254 
-0.192 ±0.135 
-0.630 ±0.194 
-0.922 ±0.179 
-0.095 ±0.164 



±0.103 
±0.209 
±0.130 
±0.268 
±0.142 
±0.203 
±0.194 
±0.172 



repeated. The differences with respect to the nominal fit are added quadratically to attain 
the systematic uncertainty. 

The unusual shape of the bins, particularly those with narrow regions, allows the possi- 
bility that a mismodeling of the invariant-mass resolution leads to an incorrect description 
of the bin-to-bin migration. In the analysis of K$K + K~ the invariant-mass resolution is 
estimated from the width of the <fi peak in data and MC. The difference is compatible with 
zero but it is conservatively assumed that the MC underestimates the resolution by up to 
one standard deviation of the measured difference. The invariant masses in the signal MC 
simulation sample are smeared by this difference and then fit to extract q and s^. This 
procedure is repeated many times and the resulting distributions of Cj and Sj returned from 
the fits are used to determine the systematic uncertainty on these parameters due to reso- 
lution. The uncertainty related to the invariant mass resolution in D° — > Kgii + 7T~ events 
is evaluated in similar fashion with the momentum of the charged tracks smeared by the 
CLEO-c resolution as in the previous analysis [7j. 

The mode-to-mode normalization in the fitter is performed either using the mea- 
sured branching fractions of the modes combined with N D0 ^o, or the measured ST 
yields. For the KgK + K~ analysis the former method is used apart from normalizing the 
Kg L K + K~ vs. Kg7i + 7T^ yields, where the pseudo-flavor ST yields are used. In the Kg7i + 7i^ 
analysis the ST yields are used apart from the two final states that contain unreconstructed 
particles in the tag candidate (K~e + v and K^ir ). In addition, previous investigations [33J 
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TABLE XVIII: Statistical and systematic uncertainties on a and Sj determined for the M = 3 
equal A6d binning of D° — > K^K + K~ data. 



1 n pnvf o i Tin r 

UllCCl idiiiLy 




Cl 




C2 




c 3 




■Si 




S2 




s 3 


(Pseudo-)flavor statistics 


n 
u 


nn^ 
uuo 


u 


UU 1 


A 
U 


UOD 


A 
U 


uro 


A 

U 


Ulo 


A 
U 


(ion 
Uoa 


Momentum resolution 


n 
u 


uuz 


u 


nn/i 

UU4 


A 
U 


U1Z 


A 
U 


m s 
Ulo 


A 

U 


UZD 


U 


UoZ 


Mode-to-mode normalization 


n 
u 


UU4 


n 
u 


nns 

UUo 


A 
U 


Ul / 


A 
U 


nni 

UUl 


A 

u 


m n 

U1U 


A 

U 


nn/i 

UU41 


Multiple-candidate selection 


n 
u 


nn/i 

UU41 


n 
u 


(1(15 
UUo 


A 
U 


U10 


A 
U 


nn/i 

UU4 


A 
U 


nns 

UUo 


A 

U 


nno 
uuz 


1 y V. k .) LUllcLLlUll 


n 


001 





om 


n 


003 


n 


009 


o 


005 


n 


009 







006 





011 





036 





132 





063 





135 


Fitter assumptions 





008 





001 





013 





013 





003 





043 


I\/T I , at" ti 1" 1 1: - !" l pc tmt* npfpTminiTio' T 1 





005 





007 





057 





024 





051 





048 


Parameterization of non-X£ final state background 





001 





001 





006 





000 





008 





003 


Parameterization of final state background 





034 





020 





061 





038 





015 





071 


Background Dalitz space distribution 





006 





015 





062 





005 





029 





022 


Assumed background B 





004 





014 





032 





001 





007 





009 


Total systematic 

Statistical plus K^K + K~ model 


0.038 
0.063 


0.034 
0.092 


0.131 
0.329 


0.142 
0.222 


0.094 
0.234 


0.175 
0.432 


K[ j K + K~ model alone 


0.000 


0.000 


0.136 


0.007 


0.000 


0.039 


Total 


0.073 


0.098 


0.354 


0.264 


0.252 


0.466 



have determined that there are small differences in the particle-reconstruction efficiency be- 
tween data and the MC. These differences lead to systematic uncertainties, and in some 
cases corrections, for final states in which K§ L h + h~ is normalized using the branching frac- 
tion and the A^^o; when using ST yields to normalize, these corrections and uncertainties 
cancel. The reconstruction efficiency of each final-state containing a 7r° or an rj meson must 
be corrected. The DT efficiency has an uncertainty for each final-state particle reconstructed 
due to the MC modeling. These uncertainties are: 0.3% per tt + , 0.6% per K + , 1.3% per n°, 
4% per rj, 0.9% per Kg, and 3% per K° L candidates. If a final state contains two particles of 
the same type the uncertainty on each identical particle is treated as fully correlated. The 
uncertainties on the ST yields, branching fractions, N D0 ^o, and particle reconstruction effi- 
ciency are used to estimate the systematic uncertainty from the mode-to-mode normalization 
method. 

The method used to select the best candidate in an event when there are multiple candi- 
dates can introduce a bias. To estimate the systematic uncertainty related to the multiple- 
candidate selection the simulation is used to determine how often the wrong candidate is 
selected. This information is used to derive corrections to the yields. The difference be- 
tween the Ci and Sj parameters fit with and without this correction applied is taken as the 
systematic uncertainty due to the multiple-candidate selection. 

The use of D° — > K~tt + , D° — > K~7t + tt° and D° — > K~tt + 7t + tt~ pseudo-flavor tags has 
a systematic uncertainty from the values of the three amplitude ratios td, the three strong- 
phase differences S d , and the two coherence parameters required to estimate the corrections 
from DCS contamination. The uncertainty is estimated by performing the fit to q and Sj 
for each parameter shifted by plus or minus the uncertainties given in Table |V} The largest 
change in each Cj and Sj, from either the positive or negative shift in the parameter, is 
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TABLE XIX: Statistical and systematic uncertainties on a and Si determined for the equal binning 
of D° -»• K° s 7r + 7r~ data. 



Uncertainty 


Cl 




C2 


C3 




c 4 


C5 


C6 


C7 


C8 


(Pseudo-)navor statistics 


0.005 





010 


0.009 





012 


0.005 


0.013 


0.013 


0.010 


Momentum resolution 


0.007 





013 


0.016 





022 


0.007 


0.021 


0.021 


0.016 


Mode-to-mode normalization 


n new 
0.00 / 





m n 
010 


0.015 





018 


n nno 
U.UU8 


n m A 
0.U14 


n no A 
U.U24 


n m 9 
0.U13 


TV 1 j_ * 1 1 * J J. 1 i ■ 

Multiple-candidate selection 


0.014 





014 


0.024 





022 


0.008 


0.014 


0.032 


0.019 


DCS correction 


n nm 
U.UU1 


U 


nno 
UUz 


U.UUl 


n 
u 


UUz 


n nno 
U.UUz 


n nn/i 
U.UU4 


n nnQ 
U.UUo 


n nnQ 
U.UUo 


LJcXIIKiZj JJ1UL clL.-L.cJJ LdllL-L. 


U.UU^i 


n 

u 


nn^i 

UUtJ 


D DDQ 


n 
u 


uuo 




D DDQ 


n 0.1 1 


n nnfi 

U.UUU 


Tag-side background 


0.024 





032 


0.049 





059 


0.027 


0.046 


0.079 


0.046 


Kgir + ir~ signal-side background 


0.014 





020 


0.028 





034 


0.016 


0.025 


0.049 


0.026 


K® 7r + 7r~ signal-side background 


0.017 





035 


0.032 





047 


0.017 


0.022 


0.046 


0.032 


Continuum background 


0.020 





026 


0.031 





038 


0.017 


0.029 


0.049 


0.031 


Total systematic 


0.042 





063 


0.080 





098 


0.042 


0.072 


0.124 


0.075 


Statistical plus K\ tt + tt^ model 


0.036 





068 


0.088 





119 


0.045 


0.102 


0.105 


0.069 


i^r7r + 7r _ model alone 


0.013 





018 


0.039 





068 


0.024 


0.040 


0.068 


0.034 


Total 

1UI Oil 


056 


n 


0Q3 


119 


n 


1 54 


062 


1 25 


n 163 


1 09 


Uncertainty 


Si 










s 4 






S7 


S8 


(rseudo-Jilavor statistics 


0.008 





014 


0.022 





021 


0.012 


0.019 


0.042 


0.017 


Momentum resolution 


0.021 





037 


0.030 





041 


0.019 


0.041 


0.039 


0.030 


Mode-to-mode normalization 


n nm 
U.UUl 


U 


nnn 

uuu 


U.UUZ 


n 
U 


nno 
UUz 


n nm 
U.UUl 


n nno 
U.UUz 


n nm 
U.UUl 


n nm 
U.UUl 


Multiple-candidate selection 


0.016 





007 


0.008 





010 


0.004 


0.007 


0.010 


0.007 


DCS correction 


0.004 





003 


0.011 





004 


0.001 


0.007 


0.016 


0.008 


Dalitz plot acceptance 


0.006 





004 


0.005 





005 


0.005 


0.004 


0.007 


0.004 


Tag-side background 


0.004 





001 


0.007 





005 


0.002 


0.003 


0.004 


0.005 


Kgir + -iT~ signal-side background 


0.002 





005 


0.005 





005 


0.005 


0.008 


0.003 


0.007 


K\i: Jr ix~ signal-side background 


0.031 





051 


0.016 





054 


0.010 


0.035 


0.074 


0.024 


Continuum background 


0.003 





004 


0.006 





006 


0.001 


0.002 


0.006 


0.004 



Total systematic 0.043 0.066 0.044 0.072 0.026 0.059 0.096 0.045 

Statistical plus ^vr+vr^ model 0.098 0.182 0.086 0.202 0.131 0.197 0.131 0.150 



KItt+tt- model alone | 0.037 0.038 0.000 0.000 0.030 0.006 0.000 0.025 

Total 0.106 0.193 0.097 0.214 0.133 0.206 0.162 0.157 



taken as the symmetric systematic uncertainty on the parameter. The total uncertainty is 
the sum in quadrature of the individual parameters. The uncertainties on the strong-phase 
differences and coherence parameters dominate the total uncertainty from this source. 

B. KgK+K "specific systematic uncertainties 

There are specific systematic uncertainties related to the K®K + K~ analysis: the strong- 
phase parameters of the Kgir + 1T~ events in the analysis, assumptions in the fitting procedure, 
the determination of the migration matrix, and the background assumptions. The use of 
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Kg L K + K~ vs. K s 7t + tt~ and K S K + K~ vs. K1tt + tt~ events introduces a dependence on the 

values of cp and sP for D° — > K® s tt + 7t~ derived from the equal A5d binning reported in 

this paper. The nominal fit has the values of dp and sP fixed to the central values. To 
determine the systematic uncertainty the input values are smeared by their uncertainties 
and the fit repeated. This is done many times and the width of the distribution of Cj 
and Si parameters returned for the particular K S K + K~ binning is taken as the systematic 
uncertainty from this source. 

The K S K + K~ fit was tested on samples of signal MC data. It was found that there 
were small but statistically significant biases in the central values of Cj and returned 
by the fitter. These biases are likely to be consequences of the assumptions used in the 
fit such as the finite granularity of the Dalitz plot bit-map and DT branching fractions 
not being corrected for quantum correlations. The whole bias is conservatively taken as 
an additional source of systematic and is found to be significantly smaller than both the 
statistical uncertainty and the dominant sources of systematic uncertainty. (No such bias is 
observed in the validation of the K s n + 7c~ fitter [7j; therefore, no uncertainty is attributed 
to fitter bias in that analysis.) 

The elements of the migration matrices, U, are determined from MC samples of signal 
events. The resulting statistical uncertainty on the elements due to the finite sample size 
introduces a systematic error on the and s^ parameters. The uncertainty is determined 
by smearing the elements by their statistical error accounting for correlations. Then the fit 
is repeated with the new migration matrices. The procedure is repeated many times and 
the width of the resulting distribution of and s^ is taken as the systematic uncertainty 
from this source. 

The background parametrization contains several sources of systematic uncertainties. For 
final states without a K\ meson the limited statistics of the generic MC sample, which is 
used to determine the peaking backgrounds, leads to an uncertainty in the background 
parameterization. For K\ there are also uncertainties arising from the ratios of signal-to- 
background events in the signal region, as well as the high and low missing-mass sidebands, 
which are obtained from the simulation. This fact, combined with the more significant 
backgrounds in the modes containing a K\ meson, leads to a larger uncertainty due to MC 
parameterizations of the background for final states that include a K\ meson. We vary these 
parameters by their uncertainties, repeat the fit then compare the result to the nominal fit 
to determine the uncertainty from this source. 

The assumptions about the distribution of the background events over Dalitz space are 
also varied to determine a systematic uncertainty. In the nominal fits the distribution 
for each background type is either assumed to be uniform or to follow the K S K + K~ signal 
distribution. We assign a systematic uncertainty to these assumptions by randomly choosing 
the fraction of the background that is uniformly distributed for each source of background; 
the remainder of that background component is then assumed to follow the distribution of 
D° — > K S K + K~ events. The fraction is assumed to have equal probability of taking any 
value between and 1. Then the fit is performed with these assumptions about the Dalitz 
plot distribution. This procedure is repeated many times and the resulting distributions of q 
and Si are used to determine the systematic uncertainty related to the assumed distribution 
of background over the Dalitz space. 

The final source of uncertainty is related to the branching fractions, £>, assumed when 
generating the generic MC samples. For modes that contribute background we vary the 
branching fraction by the uncertainty reported in Ref. [23] and repeat the fit to determine 
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the systematic uncertainty on cf^ and sf\ 

C. Uncertainties related to the acceptance and backgrounds in the K s 7v + n~ anal- 
ysis 

The sources of uncertainty related to the acceptance and background in the K s 7t + tt~ 
analysis are identical to those considered in Ref. [TJ. In addition, the same evaluation 
procedures are adopted, which we briefly outline here. 

Any difference in the relative efficiency over the Dalitz space can bias the results. In the 
nominal fit the Dalitz plot acceptance is taken from simulation. To account for any difference 
between data and simulation, the relative efficiency is smeared in each bin by 2%, which is 
the spread of efficiencies observed in simulation, and a sample of signal MC is fit with the 
new efficiencies. This procedure is repeated many times and the resulting distributions of 
Ci and Si returned from the fits are used to determine the systematic uncertainty on these 
parameters due to modeling the acceptance. (The relative variation of efficiency is an order 
of magnitude greater for D° — > K S K + K~ due to the momentum dependence of the charged 
kaon detection efficiency. Therefore, adopting a similar approach would be too conservative. 
However, the observed relative variation in efficiency from the momentum and migration 
matrix smearing for D° — » K S K + K~ is around 5%, which is greater than the fluctuations 
assumed for D° — > K s 7t + 7t~ .) 

Tag-side background yields are evaluated mode-by-mode from sidebands. The distri- 
bution over Dalitz space is assumed to follow that found in the simulation. The fits are 
repeated assuming the background events are distributed uniformly over Dalitz space. The 
difference between the fits is taken as the systematic uncertainty. 

The nominal fit ignores background in the signal K s 7r + 7c~ Dalitz space. The level of 
background estimated from K s mass sidebands is 1.9% [TJ. The background distribution is 
estimated from MC samples that include the effect of quantum correlations and the fit is 
repeated with the background subtracted. The difference between the nominal fit and the 
fit including background is taken as the systematic uncertainty due to signal-side K s it + it~ 
background. 

The nominal fit assumes the signal-side Klir + 7T~ background is distributed over Dalitz 
space following the distribution observed in the K\ mass sidebands. The fit is repeated 
assuming the background is distributed uniformly over Dalitz space and the difference in 
the parameter values taken as the systematic uncertainty. 

The continuum background is only significant in i^°7r + 7r~ events tagged by K s n° and 
K s 7c + tt~ decays. The distribution of the continuum events over Dalitz space is taken from the 
simulation. The fit is repeated with the continuum events distributed uniformly over Dalitz 
space. The difference between the fitted values of q and s, between the two distribution 
assumptions is taken as the systematic uncertainty related to the continuum background. 

D. The K^h + h~ model uncertainty 

The statistical uncertainty returned by the fit includes a contribution that is related to 



the assigned uncertainties on Aq and As« present in the \ 2 term of Eqs. (23) and (24). 
In order to isolate this contribution the fit can be repeated fixing Aq and As,, then the 
variance related to the constraint is the difference between the variances returned by the 
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fit with fixed or constrained values of Aq and As». However, due to the variations in the 
central values of Sj and q and the non-Gaussian behavior of the parameters near the edges 
of the physical region (c? + sf < 1), the covariance matrix related to the uncertainty on the 
constraint is found to give unphysical values in some cases; these are either negative diagonal 
elements or off-diagonal terms that correspond to correlations greater than one or less than 
minus one. Therefore, we present the statistical covariance matrix from the constrained fit 
and the systematic covariance matrix from all other sources of systematic uncertainty. 

To indicate that the size of the uncertainties related to Aq and As; are not dominant, 
the estimated values are given in Tables XVIII and XIX The procedure outlined above is 
used to estimate these K^lh + h~ model uncertainties. In cases where the diagonal elements 
of the covariance matrix are negative the uncertainty is set to zero. 



E. Summary of systematic uncertainties 

In both the Kgh + h~ analyses the largest sources of uncertainty is usually the background 
parameterization. In some bins there are significant contributions from the momentum res- 
olution, flavor-tag statistics, and the knowledge of strong-phase parameters. The determi- 
nations of Ci have a systematic uncertainty comparable in size to the statistical uncertainty 
whereas for measurements the statistical uncertainty dominates. Correlation matrices for 
the systematic uncertainties associated with the different binnings are given in Ref. [28] . 



VII. FINAL RESULTS, IMPACT ON 7/^3 AND CP CONTENT OF D° 



K° S K+K- 



Section 



VII A 



presents the final results for and s\'\ along with comparison of the 
measured values with the amplitude model predictions. The impact of the q and Sj results 
on the determination of 7/^3 is discussed in Sec. VII B In addition, results for the CP-odd 
fraction, J-l, of D° 

Sec, rvnci 



-> K° S K- 



K decays in the region of the <p resonance are presented in 



A. Final results for c> and S; 



The measured values of c% and Sj are shown in Figs. 13 and 14 for D° — > KgK + K 



and D° — > K^tt + tt , respectively. Also, shown are the expectations from the BABAR 2010 
model and BABAR 2008 model for D° K^K + K~ and D° K^n + 7T~, respectively. To 
test the compatibility of our results with the predictions of the models we compute 

X 2 = (P - p^°del)T V -l^p _ pmodel) ? ^ 

where P is a vector of the 2J\[ measured values of q and s,, p model is the equivalent vector 
of q and Si values predicted by the models, and V is the 2Af x 2j\f combined statistical and 



systematic covariance matrix. Table XX gives the values of the \ 2 and the corresponding 
probabilities. In the case of D° — > K®K + K~ the agreement between measured and predicted 
values of c, and Sj is good for all binnings. For D° — > Kg7i + 7i^ the compatibility of the 
model predictions with the optimal and modified-optimal binning results is very good; the 
compatibility of the models with the BABAR 2008 and Belle model equal A5o binnings is 
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FIG. 13: Measured values of Cj and Si for D° -)■ K Q S K + K' data divided into (a) M = 2, (b) M = 3, 
and (c) M = 4 equal A<5d bins. The expected values calculated from the BABAR 2010 model are 
indicated by the stars. The circle indicates the boundary of the physical region cf + s| = 1. 
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FIG. 14: Measured values of q and Sj for (a) the optimal binning, (b) the equal ASjj binning, 
(c) the modified-optimal binning, and (d) the equal A5d binning based on the Belle model. The 
expected values calculated from the BABAR 2008 model are indicated by the stars. The circle 
indicates the boundary of the physical region cf + sf = 1 . 
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reasonable. The compatibility with the predictions from the Belle model is also tested and 
is found to be reasonable. 

TABLE XX: Values of x 2 P er degree of freedom (DOF) and the corresponding probability for 
comparison of measured values of q and Si with those predicted by the BABAR 2010 and BABAR 
2008 models for D° -> K Q S K + K~ and D° -> K^tt+tt', respectively. 



Binning ^/DOF Probability 



D° - 


> K° S K + K- 




M = 2 equal A5d 


1.7/4 


79% 


M = 3 equal ASjj 


1.4/6 


96% 


M = 4 equal A5d 


2.2/8 


98% 


D° - 






Optimal 


15.5/16 


49% 


BABAR 2008 equal A5 D 25.3/16 


6.5% 


Modified optimal 


13.8/16 


61% 


Belle 


26.8/16 


4.4% 



B. Impact of Cj and Sj measurements on determining 7/03 

Our measurements of q and Sj have two consequences on the determination of 7/03 from 
a model-independent analysis: 

• the uncertainties on q and s« will result in a systematic uncertainty on the determi- 
nation of 7/«/>3, and 

• the choice of binning affects the statistical precision of a 7/03 measurement. 

Therefore, we investigate the impact on future 7/03 measurements using a simplified MC 
simulation of 5-decay data on which we perform the model-independent determination of 
7/^3 0- The number of B~ decays in a given bin of Dalitz space is dependent on q, Sj, 
r B, 5b, 7/03, and normalization parameters. The input values used for Tb, 5b, and 7/03 
are 0.1, 130°, and 60°, respectively, which are consistent with current measurements [21 E]. 
The values of are those predicted by the respective models. The yield in each bin is 

generated randomly according to the input values, and a y 2 fit between the observed and 
expected events in each bin is performed to extract the value of y/03- 

To determine the systematic uncertainty in future y/03 measurements due to the uncer- 
tainty in the measured strong phases we generate samples of five million signal events using 
the measured values of c, and Sj. The samples are then fit using the measured central values 
smeared according to their uncertainties, with the correlation between parameters taken into 
account. The large number of signal events means that the width of distribution of fitted 
7/03 gives the systematic uncertainty due to the uncertainty in Cj and Sj, as the intrinsic 
width due teethe statistical fluctuations for five million events is negligible in comparison. 

For the D° — > K$K + K~ mode, the induced uncertainty on 7/03 due to the total un- 
certainties on q and Si is evaluated from the RMS of 7 distribution returned by the fits to 
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the simulated experiments. The induced uncertainty is 3.9°, 3.2°, and 3.9° for two, three, 
and four bins, respectively. The larger uncertainty for four bins reflects the limited statistics 
available to determine the parameters in each bin. It can be noted that there is only a 
limited improvement in the statisitical sensitivity to 7/^3 by increasing the number of bins, 
since most of the sensitivity is due to the dominant resonances, K®(j) and K^a (980), which 
lie in the same region of the Dalitz plot close to the tuk+k- threshold. 

For the D° — > A"<j7r + 7r - case we find the uncertainty on 7/^3 related to the uncertainties 
on q and Si in a similar manner. The uncertainty varies significantly between the different 
binnings. For the equal A5& binning derived from the BABAR 2008 model, we find an 
induced uncertainty of 2.0°, whereas for the optimal-binning it is 3.9°. These variations 
reflect the limited statistics available with which to determine and Sj in some of the optimal 
bins; these are the same bins that lead to the reduced statistical performance in determining 
7 from the optimal binning in the presence of background. This is further emphasized by 
the expected uncertainty on 7/^3 from the modified-optimal binning results being 2.1°. The 
values of q and s» measured for the equal A5d binning derived from the Belle model lead 
to a 1.7° uncertainty on 7/^3. In all cases the error on the predicted uncertainty is less 
than 0.1° and the contribution of the systematic uncertainties on q and Sj is approximately 
1.3°. For the majority of the i^7r + 7r - binnings, these results demonstrate that, with the 
exception of the optimal binning choice, the uncertainty on 7/(^3 arising from the errors on 
the measured Cj and Sj parameters is less than that of the assigned model uncertainty in the 
unbinned analyses [T31 E] • 

For both D° ->• K° S K + K~ and D° -> K^ir + n~, small biases of 0(1°) are observed in the 
mean fitted values of 7/^3. This was also observed in Ref. [7], where the origin of the bias 
is linked to the assumed values in the fit being unphysical, c? + s? > 1, in some simulated 
experiments. If the unphysical simulated experiments are removed the bias is not eliminated 
because of the non-Gaussian nature of the truncated distributions of q and Sj. As before it 
is found that improvements in the precision of the measurements of Cj and Sj would reduce 
the bias to a negligible level. 

We also investigate how the choice of binning for D° — > Kg7T + 7T~ affects the statistical 
precision on 7/(^3. The optimal and the modified-optimal binnings have been optimized 
to minimize the statistical uncertainty on 7/(^3 in the absence of background and with 
the anticipated background at LHCb [T6], respectively. The impact of these binnings on 
B data can be assessed by comparing them to the equal A5d binning derived from the 
BABAR 2008 modeFusing simulated B data samples. We generate simplified MC samples 
of 5000 B ± K ± D° signal events which corresponds to the expected yield with 2 fb _1 
of LHCb data. The statistical uncertainties on 7/^3 are 8.0° and 7.0° for the equal A5d 
and optimal binnings, respectively. This demonstrates the advantage in using the optimized 
binning, and the improvement is consistent with the increase in the Q quantity given in 
Table |T} We also generate samples that include background events with yields according to 



the LHCb background model described in Sec. Ill B 4 We find the statistical uncertainties 



on 7/^3 of 15.0°, 15.4°, and 14.6° for the equal A5p, optimal and modified-optimal binning, 
respectively. This shows that in the presence of background the modified-optimal binning 
has a smaller statistical uncertainty than the other binnings as expected from the Q' values 
given in Table [T| 

As an additional cross check of the sensitivity to 7/^3 we determine for the different 
binnings from the measured values q and Sj. The measured values of Q' are given in 
Table XXI In all cases the values indicate good sensitivity to 7/^3 relative to the unbinned 
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TABLE XXI: Values of Q^' derived from the measured values of Cj and S{ for the different binnings. 
The uncertainty is determined by varying a and Sj within their errors accounting for correlations. 



Binning 


Q ( 





D° -»• K° S K+K- 






J\f = 2 equal A5d 


0.88^ 


-0.14 
-0.08 


M = 3 equal A<5d 


0.82^ 


-0.15 
-0.06 


Af = 4 equal A<5 D 


0.90j 


-0.21 
-0.03 


£>° -»• Kgvr+vr- 


Optimal 


0.93;! 


-0.12 
-0.02 


BABAR 2008 equal Afe 


0.8lj 


-0.05 
-0.01 


Modified optimal 


0.95j 


-0.10 
-0.04 


Belle equal A8d 


0.78^ 


-0.05 
-0.01 



method. In addition the values are in reasonable agreement with the predictions presented 
in Sec. ED 



C. CP fraction in the region of the resonance in D° 



K°K+K- 



decays 



with different criteria on Am 2 K+K 



TABLE XXII: Measured values of T- (T + ) for D° -> K%K+K~ (D° -> K\K + K~) decays 

The predicted value of AT and the average value of T- for 
and D — > K^K + K~ are also given. The hrst uncertainty is statistical and the 
second is systematic. 



D° -> KgK + K 



Criterion on 


D° -> K^K+K- 


D° -> KlK+K- 


AT 


Combined 


\&rn 2 K+K _\ 


T- 






T- 


< 0.006 GeV 2 /c 4 


1.09 ±0.09 ±0.09 


0.98 ±0.09 ±0.02 


-0.02 ±0.04 


1.03 ±0.07 ±0.04 


< 0.010 GeV 2 /c 4 


1.13 ±0.08 ±0.08 


0.98 ±0.07 ±0.02 


-0.02 ±0.04 


1.05 ±0.06 ±0.04 


< 0.014 GeV 2 /c 4 


0.73 ±0.27 ±0.04 


0.90 ±0.09 ±0.02 


-0.02 ±0.05 


0.90 ±0.10 ±0.02 


< 0.018 GeV 2 /c 4 


0.75 ±0.22 ±0.04 


0.90 ±0.08 ±0.02 


-0.03 ± 0.06 


0.90 ±0.09 ±0.02 



The CP-odd fraction of D° -> K S K + R- 



decays has been estimated using the expression 
given m Eq. ([13]. The values of Mf_ }i are measured for four different bins defined as having; 
\Am 2 K+K . 



less than 0.006, 0.010, 0.014, and 0.018 GeV 2 /c 4 . Here, Am 2 K+K . 



2 2 
m K+K-~ m <b: 



where is the nominal mass [23] . The different ranges allow the result best suited to the 
experimental resolution to be used in evaluating any systematic unc ertaint y related to the 
CP content of the D° — > K®(f) decay. The results are given in Table XXII The systematic 



uncertainty contains contributions from the migration matrix uncertainties, background 
parameterizations, and branching fraction uncertainties. The methods for determining the 
systematic uncertainties are identical to those used for c~P and sf . 

To improve the precision of the measurement we also determine the CP-even fraction, T+, 
for D° — > K®K + K~ decays in the same mass-squared intervals. The measurements are given 
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in Table XXII the significant sources of systematic uncertainty are the same as those for 
the value of T- measured from D° — > KgK + K~ decays, but are largely uncorrelated. The 
value of JF + for D° -> K L K + R- decays will be slightly different to T_ for D° -> K° S K+K+ 
decays for the same reason that q and Sj differ from c[ and s' i: as discussed in Sec. |V} 
Therefore, we correct the measured value of J- + before combining it with the measurement 
of T-. We define A J 7 as the value of T+ for D° — > K^K + K~ decays minus the value of 
T- for D° — y KgK + K~ decays predicted by the BABAR amplitude model 1141 . The value 



of AJ 77 is given in Table XXII the uncertainty on AJ 7 is determined in identical fashion to 
those on Aq and Asj, as described in Sec. |V} We then subtract the value of A J 7 from the 
measured value of T+ for D° — > K\K + K~ decays and average the result with J-_. The 
combined values of T- are given in Table XXII The systematic uncertainty on the average 
value of T- has a significant contribution from the error on AT as well as the other sources 
already discussed. 

The results are consistent with no contamination from CP-odd eigenstates in the region 
of the resonance for all m 2 K+K - intervals. The values of T- can be greater than one due 
to the background subtraction resulting in a negative yield. We calculate the lower limit 
on T- at the 90% confidence level (CL), by integrating the Gaussian distribution of the 
average value of T- within the physical region < T- < 1. The lower limits on T- at the 
90% CL are 0.89, 0.91, 0.76, and 0.77 for \Am 2 K+K . \ less than 0.006, 0.010, 0.014, and 0.018 

r2 /„4 



GeV /c 4 , respectively. 



VIII. SUMMARY 

Using 818 pb _1 of data collected by the CLEO-c detector at the -^(3770) resonance we 
have presented measurements of the amplitude-weighted averages of the cosine and sine 
of the strong phase differences beween D° and D — > Kg L h + h~ (h = 7T,K) decays in 

bins of Dalitz space, cp and sj . These results are necessary input for performing model- 
independent measurements of the CKM- angle 7/^3 [19] and can also be used in model- 
independent determinations of the charm- mixing parameters [9]. 

The results for D° — > Kg7i + tt~ are an update to those of our earlier publication [7j. The 
measurements presented here benefit from an improved analysis procedure and bin choices 
informed by an amplitude model developed by the BABAR collaboration [T2] which provides 
a better description of Dalitz space than was previously available. We have given results for 
bin choices made with an equal division in strong-phase difference ('equal AJp') and for a 
division which optimizes the foreseen precision on 7/^3 in a low background environment 
('optimal'), as expected at an e + e~ experiment, and for one ('modified optimal') which gives 
the best result under the background conditions anticipated at a hadron-collider experiment 
such as LHCb. Results have also been presented for an equal A8d binning based on an 
amplitude model devised by the Belle collaboration [13J. We estimate the uncertainty on 
7/^3 to be between 1.7° and 3.9°, depending on the binning, due to the uncertainties on the 
measured values of q and S{. In most cases, this uncertainty is smaller than that due to the 
D° — > Kg7T + 7T~ amplitude model in the most recent analyses [131 E]- 

The results for D° — > K^K + K~ are the first measurements of cp and sp for this decay. 
They have been given for equal A5d divisions of Dalitz space based on the amplitude model 
found in Ref. [H] for each half of the Dalitz plot divided into two, three and four bins. The 
uncertainty on 7/(^3 from the error on Cj and Sj parameters is between 3.2° and 3.9° depending 
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on the binnings; these uncertainties are comparable to those related to the amplitude model 
in the unbinned analysis [Hj. 

We test the compatibility of the measured values of q and Sj with those predicted by 
the amplitude models derived from flavor-tagged samples of D° — > Kgh + h~ decays. The 
agreement is reasonable in all cases, indicating that there are no large errors in the phase 
measurements provided by these models. 

In addition, we determine the CP-odd fraction, for D° — > KgK + K~ decays in the 
region of the <fr resonance. The results are given in different ranges of invariant mass squared 
about the <ft resonance. In all intervals considered T- is greater than 0.76 at 90% CL. This 
result will better constrain systematic uncertainties related to the the CP-even content when 
D° — > Kg(j) is used as a CP-odd eigenstate in an analysis. 
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IX. EPAPS ADDENDUMS 



A. Binning and model look-up tables 



The three binnings of the D° — > K^K + K~ Dalitz plot are provided as look-up tables that 
consist of grids of (m^,m 2 K+K ^), in units of GeV 2 /c 4 , and the bin number at that point. 
The granularity of the grid is 0.00179 GeV 2 /c 4 x 0.00536 GeV 2 /c 4 . The four binnings of the 



7T 



~7r Dalitz plot are provided as look-up tables that consist of grids of (m+, m\ 



in units of GeV 2 /c 4 , and the bin number at that point. The granularities of the grids for 
BABAR model and Belle model derived binnings are and 0.005 4 GeV 2 /c 4 x 0.0054 GeV 2 /c 4 
and 0.00526 GeV 2 /c 4 x 0.00526 GeV 2 /c 4 , respectively. Table |XXIII| lists the binning and 
the respective filename. 

In addition, look-up tables are also provided that contain the amplitude and phase in- 
formation for the different models over the Dalitz plot. The granularity is the same as the 
binning look-up tables. For the BABAR models 8b at each point is provided and for the 
Belle model |A5d|- Table XXIV lists the model and the respective filename. 



TABLE XXIII: Binning look-up table filenames. 



Binning 






Look-up table filename 




D° -> K° S K+K- 


Equal A5d 


M = 


2 


BinningLUT_K0SKK_BABAR2010_EqualDeltadeltaD_2bins. 


txt 


Equal A8b 


N = 


3 


BinningLUT_K0SKK_BABAR2010_EqualDeltadeltaD_3bins 


.txt 


Equal A5b 


N = 


4 


BinningLUT_K0SKK_BABAR2010_EqualDeltadeltaD_4bins 


.txt 


D° -> K^ir+ir- 



BABAR 2008 equal A5d BinningLUT_K0Spipi_BABAR2008_EqualDeltadeltaD.txt 

Optimal BinningLUT_K0Spipi_BABAR2008_0ptimal . txt 

Modified optimal BinningLUT_K0Spipi_BABAR2008_Modif iedOptimal . txt 

Belle equal A8b BinningLUT_KOSpipi_Belle_EqualDeltadeltaD.txt 



TABLE XXIV: Amplitude and phase look-up table filenames. 



Model 




Look-up table filename 


BABAR 2010 D° -> 
BABAR 2008 D° -> 
Belle D° -> Kfyr+n 


K° S K + K- model 
K^ir + ir~ model 
model 


ModelLUT_K0SKK_BABAR2010 . txt 
ModelLUT_K0Spipi_BABAR2008 . txt 
ModelLUT_KOSpipi_Belle . txt 



B. Correlation matrices 



Tables XXV 



to 



XXXVIII 



give the correlations among c( and for the different binnings. 
The correlations are split into statistical and systematic uncertainties. The correlations due 
to the uncertainty on the Ac, and As, constraints are included in the statistical correlations. 



49 



TABLE XXV: Statistical correlation coefficients (%) among cp and sp parameters for the M 
equal A6d binning of the D° ->■ K° S K + K~ Dalitz plot. 





C2 


Si 


S2 


4 


4 




4 


Cl 


-1.1 


-0.3 


-0.3 


98.8 


-l.i 


-0.3 


-0.3 


C2 




0.3 


3.0 


-1.1 


91.4 


0.3 


3.0 


Si 






7.5 


-0.4 


0.3 


99.4 


7.5 


S2 








-0.3 


2.7 


7.5 


99.5 


c ; 










-1.1 


-0.4 


-0.3 


4 












0.3 


2.7 



TABLE XXVI: Systematic correlation coefficients (%) among cp and sP parameters for the N 
equal A5d binning of the D° — > KgK + K~ Dalitz plot. 



c 2 si s 2 c\ c' 2 s[ s' 2 

ci 24.4 6.7 1.5 58.3 20.6 9.0 5.3 

c 2 1.7 0.1 22.5 60.5 5.0 5.6 

si 60.1 4.9 2.3 95.8 52.0 

s 2 5.5 4.3 58.9 87.1 

c{ 43.2 16.1 23.1 

d 2 13.6 21.9 

si 56.8 



TABLE XXVII: Statistical correlation coefficients (%) among cp and sP parameters for the N 
equal A5 D binning of the D° ->• K%K + K~ Dalitz plot. 





C2 


C3 


si 


S2 


S3 


A 


c' 2 


c' 3 


j 
s i 


j 

s 2 


4 


Cl 


0.6 


-2.9 


-1.3 


-0.6 


-0.2 


97.6 


0.6 


-2.1 


-1.3 


-0.6 


-0.1 


C2 




0.8 


1.1 


3.2 


-0.0 


0.6 


98.0 


0.4 


1.1 


3.2 


-0.0 


C3 






-0.3 


0.2 


2.5 


-2.8 


0.8 


66.4 


-0.3 


0.2 


2.4 


Si 








-2.0 


6.1 


-1.3 


1.1 


-0.1 


99.4 


-2.0 


5.8 


S2 










-3.3 


-0.5 


3.2 


0.1 


-1.9 


99.9 


-3.0 


S3 












-0.2 


-0.0 


2.2 


6.1 


-3.3 


93.7 
















0.6 


-2.0 


-1.3 


-0.5 


-0.1 


















0.4 


1.1 


3.2 


-0.0 


4 


















-0.1 


0.1 


2.2 



si -1.9 5.8 

s' 2 -3.0 
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TABLE XXVIII: Systematic correlation coefficients (%) among cp and s- parameters for the 
M = 3 equal A5 D binning of the D° ->• K° S K + K~ Dalitz plot. 



C2 


C3 


si 


S2 


S3 


4 


4 


4 


4 


4 


4 


ci 21.9 


6.4 


7.2 


-2.2 


2.4 


65.2 


19.9 


14.5 


7.5 


0.8 


2.4 


C2 


11.9 


1.2 


-15.6 


-0.5 


14.0 


49.5 


16.6 


1.2 


-14.5 


-0.8 


C3 




-0.5 


12.6 


-5.5 


7.3 


10.3 


73.5 


-0.3 


13.4 


-6.1 


Si 






38.5 


8.0 


7.4 


4.5 


-0.1 


99.4 


38.0 


7.9 


S2 








7.7 


-0.3 


-6.1 


6.0 


38.4 


95.5 


7.8 


«3 










9.9 


10.4 


2.0 


8.4 


10.6 


96.9 


ci 












41.9 


30.7 


8.5 


10.5 


9.8 


4 














35.7 


5.7 


5.8 


10.1 


4 
















0.8 


14.2 


1.3 


4 


















38.3 


8.3 


4 




















10.7 



TABLE XXIX: Statistical correlation coefficients (%) among cp and sp parameters for the M = 4 
equal A5d binning of the D° ->■ K° S K + K~ Dalitz plot. 



C2 


C3 


c 4 


si 


S2 


S3 


S4 


A 


4 


4 


4 


A 


4 


4 


4 


a 0.6 


-0.2 


-2.6 


-0.9 


0.1 


-0.0 


0.1 


98.2 


0.6 


-0.2 


-1.7 


-0.9 


0.1 


-0.0 


0.0 


C2 


0.8 


-0.1 


0.6 


1.8 


-0.7 


0.1 


0.6 


99.2 


0.8 


-0.2 


0.6 


1.8 


-0.7 


0.1 


C3 




-0.1 


0.4 


-0.3 


6.4 


0.1 


-0.2 


0.8 


98.0 


-0.1 


0.4 


-0.3 


6.4 


0.1 


C 4 






0.6 


-0.1 


0.1 


8.8 


-2.5 


-0.1 


-0.1 


63.5 


0.6 


-0.1 


0.1 


8.4 


Si 








-8.0 


8.5 


16.0 


-0.9 


0.6 


0.4 


0.5 


99.7 


-8.0 


8.5 


14.9 


S2 










-6.6 


-1.3 


0.1 


1.8 


-0.3 


-0.1 


-7.9 


99.9 


-6.6 


-1.3 


S3 












1.5 


-0.0 


-0.7 


6.2 


0.0 


8.5 


-6.6 


99.8 


1.4 


S4 














0.0 


0.1 


0.1 


6.6 


15.9 


-1.3 


1.5 


93.2 


A 
















0.6 


-0.2 


-1.7 


-0.9 


0.1 


-0.0 


0.0 




















0.8 


-0.2 


0.6 


1.8 


-0.7 


0.1 






















-0.1 


0.4 


-0.3 


6.2 


0.1 
























0.5 


-0.1 


0.0 


6.4 


^i 
























-8.0 


8.5 


14.8 


4 


























-6.6 


-1.3 


4 




























1.4 



51 



TABLE XXX: Systematic correlation coefficients (%) among cp and sp parameters for the N = 4 
equal A5d binning of the D° — > KgK + K~ Dalitz plot. 



C2 C 3 C 4 Si 


S2 


S3 


s 4 




4 


4 


4 




4 


4 




Cl 14.5 19.7 4.3 6.3 


-4.0 


5.6 


2.2 


67.7 


15.2 


20.0 


10.9 


6.7 


-3.1 


6.8 


4.3 


c 2 19.1 6.5 2.2 


-7.3 


3.7 


-1.1 


12.4 


60.1 


17.2 


9.6 


2.4 


-6.7 


4.4 


0.1 


c 3 6.9 10.9 


-13.6 


24.8 


-0.4 


18.3 


17.8 


78.5 


12.9 


11.3 


-12.6 


25.7 


1.8 


C4 1.3 


4.1 


1.3 


-1.6 


8.8 


10.6 


10.3 


80.9 


1.7 


5.0 


2.6 


-0.7 


Si 


-31.0 


71.0 


34.7 


5.4 


3.0 


9.4 


0.3 


99.7 


-30.5 


70.3 


33.2 


S2 




-27.8 


-0.1 


-1.8 


-3.3 


-9.6 


3.9 


-30.9 


99.2 


-27.2 


1.9 


«3 






1 Q Cl 

lo.y 


4.0 


Q 1 
O.l 


1 O 1 


0. / 


/U.8 


—2 1 .5 


no ^7 

y8. < 


12. o 


S 4 








2.1 


0.2 


0.4 


-6.2 


34.6 


0.1 


13.9 


96.7 


ci 










37.3 


35.1 


23.2 


7.1 


2.3 


10.0 


10.9 


c' 2 












35.4 


24.0 


4.8 


1.0 


9.1 


9.4 


c's 














23.1 


10.9 


-6.1 


23.8 


8.0 


<4 
















1.4 


6.4 


4.2 


-2.6 




















-30.2 


70.5 


33.6 






















-26.1 


3.1 


s's 






















13.9 
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TABLE XXXI: Statistical correlation coefficients (%) among cf 1 and sp parameters for the optimal 
binning of the D° — > Kgir + ir~ Dalitz plot. 
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TABLE XXXII: Systematic correlation coefficients 
optimal binning of the D° — > K^tt + tt~ Dalitz plot. 
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TABLE XXXIII: Statistical correlation coefficients (%) among cp and parameters for the equal 
A<5 D binning, according to the BABAR 2008 model, of the D° ->■ K^tt + tt^ Dalitz plot. 
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TABLE XXXIV: Systematic correlation coefficients (%) among cp and sp parameters for the 
equal A5d binning, according to the BABAR 2008 model, of the D° — > Kgir + ir~ Dalitz plot. 
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TABLE XXXV: Statistical correlation coefficients (%) among cp and sp parameters for the 
modified-optimal binning of the D° — > K^7t + tt^ Dalitz plot. 
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TABLE XXXVI: Systematic correlation coefficients (%) among cp and sp parameters for the 
modified-optimal binning of the D° — > K^7t + tt^ Dalitz plot. 
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TABLE XXXVII: Statistical correlation coefficients (%) among cp and sp parameters for the 
equal ASd binning, according to the Belle model, of the D° — > Kg-jr + ir~ Dalitz plot. 
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TABLE XXXVIII: Systematic correlation coefficients (%) among cp and sp parameters for the 
equal A5d binning, according to the Belle model, of the D° — > Kg-jr + ir~ Dalitz plot. 
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